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ABSTRACT 


This  report  deals  with  various  affects  of  the  plasma  which  forms  around  a 
re-entry  vehicle  on  antenna  characteristics^ 


A  study  of  numerical  computations- made  ory  the  IBM  704  shows  that  double-pre- 

/  \ 

cision  calculations  are  inaccurate  for  certain  ranges  of  the  plasma  parameters, 
scNthat  triple  , precision  should  be  employed.  A  linear  relation  isx  shown  to  exist 
betweei\the  change  in  input  admiitap.ce  ofa  spherical  slot  antenna  from^the  free™ 
space  veluq  and  the  change  in  retractive  index  '7la  of  the  plasma,  for  small  r\x— ( s 


/ 


as  expressed  5^  the  relation 

This  suggests  a 


\ 


AY  =  ZK(rt*-i). 

“enient  pre-flight  calibration  technique  whereby  in-flight 


This  suggests  a  convenient  pre-ilight  caliDratxon  technique 
measurements  can  bs  usea\to  deduce  plasma  properties..  \ 

mu  .  1-  ..i  .  .1 n  i.f  n. 3t _ .  .1 J  a.1.  . 


The  harmonic  series  formulation  for  the  input  admittance  is  found  to  be  im- 

\  \ 

practical  for  values  of  vehicle 'circumference  greater  than  about  one  wavelengths 

An  alternative  formulation  suitable  for  large  vehicle  sisses  1b  outlined. 


•^A  general  formulation  i3  developed  for  spherically  inhomogeneous  plasmas. 

It  is  shown  that  only  the  radial  functions  are  affected  by  the  inhomogeneity,  so 
that  results  for  the  uniform  plasma  case  can  be  extended  to  the  non-uniform  case 
by  the  substitution  of  the  appropriate  radial  functions^  -I±.~i^whowg,fIhQw«\m¥^‘‘' 
Ttasft  the  electric  and  magnetic  modes  satisfy  appreciably  different  equations  when 
the  refractive  index  gradient  in  the  _plaaaaJJ.fl.JLarge^) _ _ 


.  Investigations^of  the  effect  of'-'toe  plasma'^n  the  voltage^distri button  along 
the  slot  did  not  lead.to^a  solution  of  thls^problem.  However,  ^procedure  which 


appears  promising  is  outlined. 
^The 


effect  of  noise  generated  by  the  plasma  on  the  problem  of  signal  reception 


ii 


aboard  a  hypersonic  re-entry  vehicle  is  discussed.  It-ie-efeowrrthat  the  ef- 

it* 

fective  noise  temperature  of  the  plasma  depends  ok  the  attenuation  of  the 
plasma.  Jts^a^e6Bi^h«t®8vrthe  optimum  frequency  for  reception  usually  will  be 


significantly  lower  than  for  transmission.  An  optimization  procedure  is  de¬ 
scribed. 


iii 
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EFFECTS  OF  RE-ENTRY  PLASMA  SHEATH  OK  ANTENNA  CHARACTERISTICS 


1.  INTRODPCTICN 

This  report  summarizes  the  work  done  under  the  subcontract  on  the  theoretical 
investigations  into  the  effects  whi uh  the  plasna  sheath  generated  by  a  high  veloc¬ 
ity  vehicle  re-entering  the  earth’s  atmosphere  produces  on  the  characteristics 
of  antennas  carried  on  the  vehicle .  The  work  reported  here  constitutes  an 
extension  of  work  carried  cut  by  the  Electromagnetic  Research  Corporation  (ERC) 
and  reported  under  two  previous  subcontracts  with  the  Missile  and  Space  Vehicles 
Department  (MSVD)  [1]„ 

In  the  previous  work,  attention  was  confined  to  uniform  plasmas,  and  prin¬ 
cipally  to  strongly  ionized  plasmas.  A  strongly  ionized  plasma  acts  like  a  good 
conductor  over  a  wide  range  of  frequencies.  Consequently,  an  antenna  which  has 
been  matched  to  free  space  finds  itself  in  a  changed  environment  when  the  plasma 
forma  around  it.  This  produces  a  change  in  the  input  impedance  of  the  antenna. 

If  the  antenna  is  used  for  transmission,  this  ohange  in  input  impedance  generally 
results  in  a  reduction  in  the  amount  of  power  which  is  absorbed  from  the  source 0 
In  addition,  losses  are  introduced  by  the  attenuation  of  waves  transmitted 
through  the  plasma. 

In  a  previous  report  [2]  closed  form  expressions  were  obtained  for  the  radi¬ 
ation  properties  of  several  types  of  antennas  when  surrounded  by  a  strongly 
ionized  plasma  sheath.  Particular  attention  was  devoted  to  the  analysis  of  a 
slotted  sphere  type  of  antenna,  since  this  geometry  is  the  simplest  one  for  anal¬ 
ysis  and  yet  represents  a  fairly  good  approximation  to  the  shape  of  practical 
re-entry  vehicles.  One  of  the  tasks  under  the  present  subcontract  was  to  gener¬ 
alize  the  previous  work  to  cover i 
*  Numbers  in  brackets  refer  to  References  on  p.  go. 
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{a}  weaker  ionization  densities 

(b)  higher  radio  frequencies 

(c)  noa-syametric  geometries „ 

It  was  anticipated  that  this  would  require  a  digital  computation  program,  so 
that  the  objective  was  to  cast  the  analysis  in  a  form  suitable  for  machine  pro¬ 
gramming  o 

In  a  previous  report  under  this  subcontract  [3]  earlier  work  for  the  slotted 
sphere  antenna  was  put  in  a  form  suitable  for  programming  on  the  IBM  704  com¬ 
puter  o  Programming  and  calculations  were  carried  out  by  MSVD  for  several  fre¬ 
quencies  and  for  a  wide  range  of  plasma  ionization  densities 0  Tabulations  of 
the  numerical  results  were  delivered  to  ERG  during  the  last  month  of  the  contract 
period,,  Examination  of  these  results  led  to  the  recognition  of  several  proper¬ 
ties  characteristic  of  the  computations,,  One  of  these  was  a  particularly  simple 
type  of  relation  between  the  change  in  input  admittance  of  the  antenna  and  the 
plasma  characteristics  for  weak  ionization  densities „  Another  was  that  the  se¬ 
ries  type  of  formulation  used  was  not  suitable  for  antennas  (spheres)  large  rel¬ 
ative  to  the  wavelength,  since  the  number  of  terns  in  the  series  required  in 
order  to  obtain  adequate  convergence  turned  out  to  be  too  large o  The  results  of 
the  numerical  computational  program  and  its  limitations  are  discussed  in  detail 
in  Sec„  2  of  this  report,, 

The  properties  exhibited  by  the  computed  results  led  to  additional  theo¬ 
retical.  investigations  which  are  presented  in  Sec.  3«  An  analytical  in¬ 
vestigation  of  weak  plasma  densities  led  to  a  rather  remarkable  result „  It 
was  found  that  for  a  spherical  geometry  and  a  uniform  plasma,  the  change 
in  input  admittance  of  the  antenna  when  the  plasma  forms  is  proportional  to 
the  change  in  refractive  index  of  the  plasma  sheathe  The  constant  of  pro- 
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portionallty  is  in  the  nature  of  a  structure  constant  dependent  on  the  electrical 
dimensions  of  the  vehicle c  The  type  of  result  is  quits  general  and  is  not  limited 
to  small  vehicle  sizes  relative  to  the  wavelengtho  From  the  way  in  which  the 
linear  proportionality  was  obtained,  it  is  felt  that  a  similar  result  would  be 
obtained  with  other  regular  (separable)  geometries «  In  addition,  it  Is  suspected 
that  the  same  type  of  proportionality  probably  exists  for  practical  vehicle  geom¬ 
etries  »  This  leads  to  a  convenient  and  important  technique  for  calibrating  an 
antenna  on  the  ground  so  that  in-flight  measurements  can  he  used  to  deduce  plasma 
properties.)  The  details  are  given  in  Sec0  3» 

Sec a  3  also  contains  the  outline  of  an  alternative  formulation  of  the  input 
impedance  problem  which  should  be  suitable  for  large  vehicle  sizes®  Further 
work  along  the  lines  indicated  therein  should  form  the  basis  for  a  computational 
program  which  should  bridge  the  gap  for  high  frequencies  for  which  the  earlier 
series  formulation  was  found  to  be  inadequate® 

In  the  previous  analyses  it  was  assumed  that  the  plasma  is  of  constant 
thickness  and  uniform  density®  This  is  a  highly  idealized  assumption  since, 
in  fact,  the  plasma  sheath  around  a  re-entry  vehicle  has  large  gradients  of 
ionization  density®  The  extension  of  the  analysis  to  take  into  account  these 
gradients  is  carried  out  in  Sec®  4»  It  is  shown  that  it  is  necessary  only  to 
replace  the  radial  functions  in  the  input  admittance  by  functions  appropriate 
to  the  actual  variation  of  plasma  density.  It  is  pointed  out,  however,  that 
care  must  be  exercized  not  to  overlook  terms  which  are  usually  neglected  in 
treatments  of  propagation  in  inhomogeneous  mediums,  since  these  terms  are  not 
negligible  in  re-entry  plasma  problems. 

The  analysis  of  the  input  admittance  of  a  slot  type  antenna  involves  the 
voltage  distribution  along  the  slot.  For  a  slot  which  is  electrically  short, 
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the  voltage  distribution  is  approximately  triangular  in  shape,  both  in  free  space 
and  in  the  presence  of  the  plasma..  For  frequencies  high  enough  so  that  the  slot 
is  no  longer  electrically  short,  however,  the  assumption  of  a  triangular  voltage 
distribution  is  probably  inadequate,  and,  in  addition,  the  voltage  distribution 
may  change  when  the  plasma  forms „  Consequently,  a  method  either  for  determining 
the  voltage  distribution  or  for  calculating  the  input  admittance  without  using 
assumptions  regarding  the  voltage  distribution  would  be  desirable*  Sec.  5  of 
this  report  gives  the  results  of  some  work  carried  out  on  this  problem.  A  com¬ 
plete  solution  has  not  been  obtained,  but  a  method  of  attack  which  may  prove 
successful  is  presented. 

In  Sec.  6  the  receiving  problem  aboard  a  re-entry  vehicle  Is  considered.  In 
particular,  the  effect  of  noise  generated  in  the  plasma  upon  the  reception  problem 
is  discussed.  It  is  shown  that  the  effective  noise  temperature  of  a  receiving 
system  depends  not  only  on  the  temperature  of  the  plasma,  but  also  on  the  attenu¬ 
ation  through  the  plasma  at  the  frequency  in  question. 

A  brief  summary  of  the  work  encompassed  by  the  report,  together  with  con¬ 
clusions  and  recommendations  for  further  studies,  is  given  in  Sec0  7. 


2o  PROGRAM  OF  NUMERICAL  COMPUTATIONS 


2„1  Introduction 

In  a  previous  report  [3] *  an  earlier  solution  for  the  input  admittance  of 
a  spherical  slot  antenna  was  recast  in  a  form  suitable  for  numerical  computation.. 
This  form  was  programmed  for  the  IBM  704  computer  by  MSVD0  Computations  were 
carried  out  for  a  wide  range  of  plasma  densities  at  each  of  three  frequencies: 

14,  240,  and  3000  me.  Tabulations  of  the  results  were  furnished  to  ERG  for  study . 
The  results  of  this  computation  program  will  be  discussed  in  this  section.. 

In  [3]  it  was  shown  that  it  is  preferable  to  calculate  the  difference  in 
input  admittance  of  the  antenna  in  the  presence  of  the  plasma  (denoted  by  Y) 
from  that  in  some  reference  situation.  The  reference  situation  chosen  in  [3]  is 
that  where  the  dielectric  insulating  layer  surrounding  the  antenna  (and  over 
which  the  plasma  forms)  is  infinitely  thick  (for  which  the  input  admittance  is 
denoted  by  ) .  This  choice  of  reference  was  made  because  the  expressions  in 
the  input  admittance  simplify  considerably.  The  incremental  admittance  then 
becomes 


m&Pn*,(ofFnXsf  +  FZrXofGnXg 


] 


(2d) 


The  properties  of  the  plasma  enter  only  into  the  quantities  Fn  and  which 
depend  on  the  electrical  dimensions  of  the  antenna.  In  these  quantities  the 


properties  of  the  plasma  enter  through  its  complex  refractive  index  Yl%,  where 


n*  = 


in  which 


(2.3) 

(2.4) 


■f  being  the  frequency,  and  and  p  the  ionization  density  and  collision 
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frequency,  respectively,  in  the  plasma 0  Values  of  Y  -  then  can  be  com¬ 
puted  for  various  plasma  characteristics  by  programming  the  computations  for  a 
set  of  values  of  n,g=  If  one  of  these  values  is  n8=  1,  then  for  this  value  the 
quantity  obtained  is  ^  ,  where  Yp  is  the  input  admittance  in  free  space. 

Consequently  the  change  in  input  admittance  upon  formation  of  the  plasma  is  ob¬ 
tained  by  subtracting  this  value  from  each  of  the  other  end  values  Y~  Yc  « 

AY*  Y-Yp  =  (Y-Yx)  ~  Ofr->£).  (2o5) 

This  procedure  was  adopted  since  it  appeared  to  reduce  the  complexity  of  the 
machine  program* 

In  (2*1)  the  summation  over  a.  is  shown  as  extending  to  infinity*  In  [3], 
computation  only  up  to  J1  =  30  was  recommended,  since  pilot  desk  calculator  com¬ 
putations  indicated  that  this  was  sufficient  to  obtain  convergence  for  the  high- 
frequency  situations  of  interest  in  the  earlier  studies.  A  procedure  was  indi¬ 
cated,  however,  whereby  the  computation  of  F„  and  ,  the  quantities  which 
ultimately  control  the  convergence  of  the  series,  can  be  extended  to  larger  values 
of  n. 

In  this  section,  the  numerical  results  obtained  on  the  IBM  704  by  MSVD  will 
be  discussed,  and  four  characteristics  of  these  results  will  be  pointed  out  and 
discussed* 


2,2  IBM  704  Calculations 

2o2,l  Outline  of  Investigations 

In  [3],  a  recommended  procedure  was  given  for  the  calculation  of 
input  admittance  in  accordance  with  (2,1) »  To  check  the  programming,  initial 


*  The  programming  was  done  by  Mrs,  Ruth  Lyon  of  MSVD, 
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computations  were  made  of  the  spherical  Bessel  functions  occurring  in  Fn  and  G„* 
These  computations  involve  the  use  of  recursion  formulas  to  generate  the  func¬ 
tions  of  various  orders*  Both  ascending  and  descending  recursions  are  used,  de¬ 
pending  on  the  argument.  For  intermediate  arguments,  a  check  can  be  made  by- 
using  both  procedures *  It  was  found  that  huge  discrepancies  resulted  when  or¬ 
dinary  single-precision  arithmetic  was  used*  Pilot  desk  computations  by  ERG 
showed  that  the  recursion  method  was  very  sensitive  to  round-off  error*  Since 
the  704  carries  only  eight  significant  figures,  the  amount  of  round-off  error 
incurred  in  these  computations  using  single-precision  arithmetic  Introduced  pro¬ 
hibitively  large  errors*  Accordingly,  a  double-precision  routine  was  recommended 
a3  being  necessary,  with  the  possibility  that  triple  precision  might  be  required* 

After  the  program  was  revised  for  double-precision  computations,  check  values 
of  the  spherical  functions  agreed  to  four  or  five  significant  figures.  Conse¬ 
quently,  all  subsequent  computations  were  double  precision*  Evidence  will  be 
shown  below,  however,  that  for  certain  regions  of  plasma  properties  (i<,e*, W*} 
the  resulting  accuracy  deteriorates,  so  that  triple  precision  probably  is  needed 
in  these  regions* 

The  frequencies,  vehicle  size,  and  other  parameters  for  which  calculations 
were  carried  out  were  specified  by  Dr*  W*  C*  King  of  KSVD*  Frequencies  of  14, 

240,  and  3000  me  were  specified*  Values  of  x2  and  z  were  selected  in  a  grid  to 
cover  an  extremely  wide  range*  Subsequently  a  number  of  intermediate  points 
were  calculated,  especially  for  14  me,  for  which  the  most  complete  coverage  was 
made* 

As  stated  in  Sec*  2*1,  calculations  were  made  of  incremental  admittance  in 
accordance  with  (2*1),  where  Yz  represents  the  admittance  when  the  slotted 
sphere  is  covered  with  a  dielectric  layer  of  infinite  thickness*  Plots  of  Y-Yz 
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did  not  prove  to  be  especially  revealing,  since,  for  example,  the  imaginary 
component  reversed  sign  for  certain  values  of  When  the  incremental  ad¬ 
mittance  AY  as  determined  by  (2.5)  was  plotted,  however,  certain  revealing 
characteristics  were  found*  In  particular,  the  following  four  features  became 
evident* 

(a)  A  linear  proportionality  between  AY  and  x2»  for  constant  %  and  suffi¬ 
ciently  small  x2; 

(b)  an  interference  phenomenon  at  a  value  of  X2  near  2; 

(c)  a  degradation  of  accuracy  in  the  imaginary  part  of  AY  for  large  values 
of  x*  and  z  ; 

(d)  apparently  erroneous  values  of  AY  at  certain  "pathological"  combinations 
of  x2  and  z . 

These  properties  and  their  significance  will  be  pointed  out  in  the  discussion 
of  some  plots  of  AY  to  follow* 

2.2.2  Dependence  of  AY  on  Plasma  Properties 

Fig.  2-1  is  a  plot  of  the  real  part  of  Ay  vs.  xa  for  several  values 
of  z  for  a  frequency  of  14  me.  It  can  be  seen  that  a  linear  proportionality  ex¬ 
ists  between  x3  and  both  the  real  and  imaginary  parts  of  Ay  for  x*  below  some 
upper  value.  This  upper  value  depends  on  and  increases  with  the  value  of  2. 
Furthermore,  a  regular  dependence  on  z  also  is  evident. 

Fig.  2-2  i3  the  corresponding  plot  of  the  Imaginary  part  of  A/.  For  values 
of  x*  below  about  1,  the  same  type  of  linear  increase  with  x*  and  a  regular  de¬ 
pendence  on  z  is  evident.  For  large  x2,  however,  the  curves  become  quite  irregu¬ 
lar.  In  contrast,  the  curves  of  Fig.  2-1  remain  very  smooth  and  regular  in  this 
region.  This  feature  of  Fig.  2-2  will  be  discussed  in  Sec.  2.2.4° 
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To  show  the  dependence  on  z  in  the  linear  range  sore  clearly,  AY  was 
plotted  against  z,  for  a  constant  value  of  x®  =  0„3»  as  shown  in  Fig.  2-3*  It 
was  easy  to  see  that  a  fit  to  the  imaginary  part  of  AY  was  given  by  a  relation 


of  the  form 

(2.6) 

and  the  real  part  by 

(2.7) 

where  C  is  a  constant. 

The  solid  curves  in  Fig. 

2-3  are  plots  of  <fl  and  <£  in 

accordance  with  these  relations »  Since  Figs*  2-1  and  2-2  show  that  C  is  propor¬ 
tional  to  x2  for  x2  not  too  large,  and  in  view  of  (2*2),  the  good  agreement  of 
the  plotted  points  with  those  curve's  was  very  suggestive  of  the  relation 

xK'K-l),  (2.8) 

where  K‘  is  a  Constant*  For  14  me.  Fig.  2-3  giveB  K'~  1.4  x  10_2o 

Thus  we  arrive  at  the  relation  (2*8)  for  Af  for  sufficiently  small  ~  | . 
Since  this  region  of  7?f-l  is  synonymous  with  n^asl,  (2.6)  may  be  written  as 
AY  a  J.  2.XXnz  -  0  -J.K  Cnz-f)  =  x/<V,  (2.9) 

where  K-  2 Ket  and  V  =  ■n2-! .  It  will  bt  shown  analytically  in  Sec.  3,  in  fact, 
that  AY  is  given  by  (2.9)  for  small  V  for  any  frequency  and  vehicle  dimensions, 
K  being  dependent  on  these  parameters. 

A  noteworthy  departure  from  (2.8)  vaB  found  for  the  values  of  AY  corres¬ 
ponding  to  z  =  0  in  Fig.  2-1.  This  is  classed  as  a  "pathological19  condition, 
which  will  be  discussed  further  in  Sec.  2.2.5. 

2.2.3  Interference  Phenomenon 

In  Figs.  2-1  and  2-2  a  sharp  increase  in  AY  takes  place  near 
X2  =  2  for  small  values  of  Z.  This  produces  a  hump  in  the  curve,  which  is  more 
pronounced  the  smaller  the  value  of  Z.  This  region  is  plotted  in  greater  detail 
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in  Figs •  2-4  and  2-5=  The  peak  of  the  hump  takes  place  at  a  value  of  xs  which 
is  somewhat  dependent  on  the  value  of  z> 

This  phenomenon  is  suggestive  of  interference  between  the  waves  reflected 
from  the  inner  and  outer  boundaries  of  the  plasma  sheathe  Although  for  x*  >  1 
the  plasma  represents  a  barrier  or  attenuating  layer,  its  electrical  thickness 
is  small  enough  so  that  the  attenuation  (for  vanishingly  small  Z,  for  example) 
of  a  round-trip  passage  through  the  sheath  is  quite  low.  A  notable  exception 
to  the  general  trend  of  the  other  curves  is  exhibited  by  the  curve  for  z  =  0. 

As  mentioned  above,  this  curve  appears  to  fall  in  the  category  of  "pathological" 
values  to  be  discussed  below  in  Sec.  2.2.$. 

The  situation  is  analyzed  qualitatively  in  Sec.  3.3  hy  means  of  a  plane- 
wave  analysis. 

2.2.4  Reduced  Accuracy  of  Imaginary  Part  of  AV  for  Large  z  and  xa 
Fig.  2-6  shows  an  extension  of  Fig.  2-3  to  larger  values  of  z. 

The  solid  lines  are  the  plots  of  (2.6)  and  (2.7),  while  the  computed  values  are 
shown  by  the  encircled  points.  Clearly  the  points  for  depart  increas¬ 

ingly  from  the  curve  as  z  increases.  It  will  new  be  shown  that  this  should  be 
interpreted  as  an  error  in  the  computed  values. 

For  JS  large,  (2.2)  approaches 

nf-l  ^  £(7l2“l)  =  2v>  ='*22 

Thus  for  z»xz  »  V  becomes  very  small,  decreasing  ass  increases.  Now,  as  men¬ 
tioned  earlier,  the  analysis  in  Sec.  3  will  show  that 

A y  =  +  0(v»). 

Thus  the  computed  points  should  approach  the  curves  more  closely  as  z  increases. 
Since  the  opposite  trend  takes  place,  the  computations  become  suspect. 
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It  should  be  pointed  out,  in  view  of  the  phenomenon  to  be  discussed  in  Sec. 
2.2.5,  that  for  z»  x*»  the  real  part  of  71* -I  becomes  negligibly  small  relative 
to  the  imaginary  part.  In  order  words,  - 1  approaches  a  pure  Imaginary  quan¬ 
tity,  while  n£  approaches  unity. 

It  has  been  pointed  out  already  that  the  curves  of  A  (/si/)  in  Fig.  2-2  be¬ 
come  irregular  for  large  <z.  It  is  believed  that  the  curves  are  incorrect  in 
these  regions,  and  to  indicate  this  they  have  been  dotted.  As  further  evidence 
of  this  inaccuracy,  the  values  to  be  expected  on  the  basis  of  (2.6)  are  shown 
by  the  dashed  lines,  which  have  been  drawn  up  to  the  approximate  limit  of  X2 
over  which  this,  relation  is  expected  to  hold  (xVfe»0.3) .  As  z  Increases,  the 
computed  values  depart  more  and  more  widely  from  these  expected  values.  This 
lends  support  to  the  conclusion  that  the  dotted  regions  represent  erroneous  val¬ 
ues  . 

The  reason  for  this  anomalous  behavior  of  the  computed  values  of  is 

not  known.  The  apparent  discontinuities  in  certain  regions  may  be  due  to  a 
change  in  formulation  in  the  program  which  we  understand  was  introduced  for 
certain  reasons.  Why  this  should  affect  c&(AY)  and  not  <32.(AY)  we  do  not  know. 
The  anomalous  behavior  suggests  that  double-precision  computations  on  the  704 
are  not  adequate  for  this  problem. 

2.2.5  Pathological  Values 

It  may  have  been  noticed  already  in  Jig.  2-1  and  2-2  that  several 
of  the  plotted  points  are  "out  of  line"  with  respect  to  the  smooth  curve  through 
the  remaining  points.  Each  of  these  points  is  encircled  and  joined  by  a  verti¬ 
cal  arrow  to  the  curve  on  which  it  belongs.  Curiously,  every  one  of  these 
"pathological"  points  occurs  when  X*  —  I  +  »*. 
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Furthermore,  it  was  found  that  when  a  computation  was  made  for  this  combination 
of  x*  and  a*,  a  “pathological"1  point  always  developed .  \ 

It  was  suspected  that  same  error  had  been  made  in  the  programing  which  re- 
suited  in  this  pathological,  behavior .  However,  a  thorough  check  by  Mrs.  Lyon 
of  MSVD  showed  that  no  such  anomaly  in  the  program  existed*  It  will  now  be  shown 
that  round-off  error  is  the  suspected  cause  of  these  pathological  values. 

From  (2.2),  tia  becomes  for  these  values 

vtj  »  -in, 

which  is  pure  imaginary.  As  noted  in  the  preceding  section,  the  computed  values 
of  cSr(AY)  become  increasingly  in  error  ae  yjJ  approaohen  unity  for  small  71*— J. 
The  behavior  just  discussed  indicates  that  this  same  type  of  error  occurs  for 
imaginary  values  of  Yl* . 

Another  region  of  anomalous  behavior  is  the  curve  for  z=  0  in  Fig.  2-1. 

For  small  Xa,  the  values  should  conform  to  (2.8),  as  shown  by  the  results  of  the 
analysis  to  be  given  in  Sec.  3.  However,  for  small  Xs-,  the  real  part  of  the 
computed  value  of  AY  becomes  negative  for  an  =  0,  in  contradiction  to  the  positive 
sign  required  by  (2.8)  and  the  theoretical  analysis. 

From  the  above  behavior,  it  is  suspected  that  the  recursion  calculation  of 
the  spherical  Bessel  functions  becomes  especially  sensitive  to  round-off  errors 
for  Imaginary  and  real  values  of  'Hz.  Consequently  it  appears  that,  to  achieve 
accuracy  for  all  values  of  argument,  a  triple-precision  routine  would  be  required 
on  the  IBM  704.  The  anonalous  behavior  of  d£(AY)  discussed  in  Sec.  2.2.4  also 
suggests  that  triple  precision  is  necessary  for  this  problem. 

2.2.6  K  rj  f 

It  was  shown  in  Sso.  2o2„2  that,  for  sms'll  V ■  7i2~ !  the  values 
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of  AY  obeyed  the  relation 


AY  «  jLKi/. 

The  dependence  of  the  constant  K  on  frequency  will  now  be  examined. 

From  plots  of  the  computed  values  at  the  three  frequencies,  the  values  of 
|/Cl  in  Table  I  were  determined. 


TABLE  I 

—  M  _ 

14  7.0*10-3 

240  6.5*10“3 

3000  7.2°10“2 

The  trend  exhibited  by  this  table  is  for  to  increase  with  f .  However, 
there  is  practically  no  change  between  14  and  240  me.  In  addition,  there  were 
reversals  in  the  sign  of  for  small  values  of  e  at  240  mo,  while  at 

3000  me  the  signs  were  all  negative. 

As  a  result  of  this  unexpected  behavior  of  A",  an  examination  was  made  of 
the  values  of  F„  and  which  had  been  calculated  as  a  subroutine  and  printed 
out  so  as  to  be  available  for  study.  It  was  found  that  convergence  had  not  set 
in  at  a  value  of  n  =  29  for  f  =  240  and  3000  me.  It  appeared,  however,  that  a 
value  of  n  =  29  was  just  about  adequate  for  f  =  14  me. 

From  a  study  of  the  convergence  properties  of  (2.1)  as  exhibited  by  the 
calculated  values  of  Fn  and  it  is  inferred  that  calculations  would  have  to 
b^  carried  out  to  about  >1  =  500  for  f  =  240  me  and  to  about  rt  =  7000  for  f  =  3000 
me.  Consequently  the  use  of  (2.1)  for  these  frequencies  be cooes  impractical. 
This  is  the  condition  where  the  sphere  is  large  relative  to  the  wavelength,  under 
which  condition  the  harmonic  series  analysis  becomes  impractical.  A  method  of 
dealing  with  this  situation  is  outlined  in  Sec.  3.3. 


19 


2,3  Summary 


From  a  study  of  the  computations  made  by  MSVD,  the  following  four 
properties  of  the  results  have  been  recognizedi 

(a)  a  linear  proportionality  between  the  change  in  input  admit¬ 
tance  upon  formation  of  the  plasma,  for  weak  plasmas,  ae  expressed  in  (2,8), 

(b)  an  interference  phenomenon  for  small  z  due  to  the  waves  re¬ 
flected  at  the  inner  boundary  of  the  plasma, 

(c)  impaired  computational  accuracy  for  plasmas  characterized  by 
imaginary  or  real  values  of  the  square  of  the  refractive  index, 

(d)  unsatisfactory  convergence  properties  of  the  series  formula¬ 
tion  (2,1)  for  vehicle  sizes  large  relative  to  the  wavelength. 

In  Sec,  3,  theoretical  investigations  are  made  of  (a),  (b)  and  (d).  The 
linear  proportionality  in  (a)  is  demonstrated  quantitatively,  while  a  qualitative 
explanation  of  (b)  is  given  by  means  of  a  plane-wave  analysis,  A  procedure 
whereby  the  limitation  (d)  can  he  overcome  is  outlined.  It  appears  that  limita¬ 
tion  (c)  requires  the  use  of  a  triple-precision  program  in  order  to  obtain  satis¬ 
factory  accuracy  with  the  704  for  all  values  of  refractive  index  of  the  plasma. 


3»  THEORETICAL  EXTENSIONS  FOR  A  UNIFORM  PLASMA 


In  this  Section  several  extensions  to  the  analytical  results  obtained  in 
previous  reports  which  were  obtained  under  the  present  contract  will  be  reported. 
These  extensions  include 

(a)  the  analysis  of  weak  plasmas,  including  a  discussion  of  non-spher- 
ical  and  non-symmetrical  geometries; 

(b)  a  discussion  of  the  Interference  phenomenon  pointed  out  in  Sec, 

2,2,3; 

(c)  a  formulation  applicable  for  the  condition  where  the  sphere  is 
large  relative  to  the  wavelength. 


3«1  Weak  Plasmas 

In  Sec,  2,2,2  it  was  shown  that  the  computations  indicated  that,  for 
small  AY  obeyed  the  relation 

AY- (3,1) 
where  K  depends  on  the  frequency,  but  not  on  the  plasma  properties.  In  this 
Section,  the  relation  (3«l)  will  be  derived  analytically  for  the  case  of  a  spher¬ 
ical  geometry.  Its  applicability  to  other  geometries  then  will  be  discussed. 

The  method  of  establishing  this  result  will  be  to  show  that  for  each  value 
of  the  summation  index  n  in  (2,1)  the  corresponding  term  in  AY  y  which  will  be 
designated  by  £>Yh  >  contains  y  as  a  factor. 

As  shown  in  [3],  Fn  and  <qn  are  given,  respectively,  by 


Fn  = 


rB± 


(3,2) 


(3,3) 


in  which  /^nH1  is  the  reflection  coefficient  of  the  n**  magnetic  mode  at  the 
boundary  r  -  by  and  is  the  corresponding  reflection  coefficient  of  the 
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electric  mode*  Expressions  for  and  were  given  in  [3],  and  will  be 
rewritten  below  in  a  revised  notationo 

Rn  and  appear  in  the  expression  for  Y~  Yx  given  in  (2.1).  The  cor¬ 
responding  quantities  which  appear  in  )£-)£  may  be  denoted  by  and  , 

and  the  corresponding  reflection  coefficients  by  and  It  will  be 

shown  that 

Rn**  '  +  0  (Y%  (3.4) 

and,  similarly,  that 

Rf*  -  •  -X  *  0(y*),  (3.5) 


where  and  do  not  involve  the  plasma  properties.  These  lead  to  the  cor¬ 
responding  relations 


Fn~Rn^ 

(3.6) 

(3.7) 

so  that 

AYn  *  ji^Kn  y  OCY-). 

(3.8) 

Consequently 

AY  *  ILYn  ~ 

(3.9) 

which  is  the  result  to  be  demonstrated. 


In  older  to  make  the  analysis  least  cumbersome,  it  was  found  advantageous 
to  introduce  a  simplified  notation.  The  functions  involved  are  the  spherical 
Bessel,  Neumann  and  Hankel  functions,  and  the  arguments  of  these  functions  that 
will  be  encountered  are  k,b ,  A*  4  ,  k9b  ,  kzC  ,  and  ktC  .  These  arguments  will 
be  denoted  by  the  subscripts  1,  2,  3,  4,  5  respectively.  Since  and  can 
be  written  entirely  in  terms  of  the  spherical  Bessel  (/*),  Neumann  (7T*)  and 
Hankel  (A)f)  functions,  the  subscript  n  will  be  dropped.  Consequently,  for  ex¬ 
ample,  JnOubi  is  denoted  by  j,  ,  etc.  Thus  the  following  scheme  is  used: 
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JntAlb)  —*■  j\ 

n„  CA,b)  — ►  7t, 

hn'Ooh)-*  h, 

— ►  Jg 

Tin  (kabi  -+■  71* 

hPlktb)—*  A, 

JWM  -*■  b 

Tin  U(,b)->  7lt 

A, 

Tin  (A 71+ 

CW)-*  h4 

Tln(A^)  Tig 

fin^(kgt) bf 

The  detail*  will  be  carried  through  for  #nM,  and  then  the  corresponding  re- 

B4 

suit  written  down  for  Rn  by  inspection « 

Fran  [3]  we  find  that  may  be  written  in  the  detei-mlnantal  form 


R-1A*  0 

fi-jlCa.  C 

-  a  —  ;  _ 

Am  0 
C*  P 

N 

ie  *1 

|/»  p 

1?  ? 

I 


where 


Irt* 

4v< 

lii 

*j/| 

H 

bJi 

I*  cu\ 

|Ja 

K 

brf 

1  h,  <*A»| 

■it 

c*- 

P  »  1  I 

F  ■* 

J 

r*»a 

|itf  c*n+| 

-n. 

hsJ^j 

The  corresponding  expression  for  the  free-space  condition  ia 

,  k/  <n 


i- 


la'  e' 


where 


The  essential  step  is  to  develop  the  functions  with  arguments  A»b  and  4*C 
(subscripts  2  4»  respectively)  in  a  Taylor 8 s  series  expansion  in  terms  of  the 

corresponding  functions  of  arguments  kgb  and  k%C  (subscripts  3  and  5»  respec¬ 
tively)  o  Since 
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k%  *  ntk3  *  kt  +  (/TjflMa  =  /30+>'), 

we  write,  for  example 

j*  -*■  +  v^;M+  —  y3  4-  yhA'  +  0(>>% 

*  JS +  ^PSp*  “  *]*  +  °fV*>  *  *'  ’♦  +  OCV*')» 


where 


-  rtfn+O  , 


With  transformations  of  this  type,  we  obtain 
Ag9A^  +  yb3Ai  +  Oiv  x) 

Ci'Cf+v^Cs  4-  0<v9 
a  *  &  +  VCS&3  ♦  Ofv1) 

P  *  dF  +  wCsPs  +  O^V*) 

£■  *  £r+  +  CXv*) 

F  -F'V  vi*/*  +  OCv*) 

Then  the  difference  between  and  becomes  a  imply 

ftfL  -  *  -*v%H  +  06**), 

where 


003* 

<z  * 
€  * 


*1*7' 

‘■P7 

U*  *3 

.  4  k  ^ 

Kf  ^ 

«  1C* 

k  a. 
k  ** 

+  £  It  y 

AS  s' 

of  ^ 

24 


A  similar  analysis  may  be  performed  for  "the  electrio  modes .  The  result  may 
be  written  down  at  once  by  inspection,  since  only  the  values  of  ,  CL  and 


t  are  changed! 

*„el  '  ^  *  -vcV-7^e  +  OfV‘). 

Finally,  there  remains  the  incorporation  of  these  results  into  the  corre¬ 
sponding  differences  Fn  -F* bxA  Qn  -  6„p,  which  are  the  quantities  which  appear  in 


the  Incremental  admittance 


AY  *  /-  YF  •  f)A Yn. 

For  these  we  have,  on  denoting  functions  of  argument  kta.  by  the  subscript  6, 

r  r.  *r _ aT  ,, 


+  CXV*) 


+  cyv>*) 


'n  T  /  y  *n  ;  -  - 1 - 

LAo -,£>(1  -Jo ( I  -viA*')]  U7lofi(A«\>^%]£>'ru  UA" 


(•not-  A/<,AH,,)t 


+  0(VV). 


Similarly, 

G>«-  SrT’ 


(**  +  i#A**)‘ 


OCv1), 


Since  V  is  a  factor  in  —  endin'" for  each  value  of  n,  it  is  a  factor 
in  AY,  so  that  we  have,  in  general,  for  any  electrical  dimensions  of  the  spher¬ 
ical  antenna 


AY-  M  +  0<v»), 

which  is  (3.9),  which  we  set  out  to  derive.  Therefore,  for  sufficiently  small 
V,  A Y  is  proportional  to  V, 


3.1.1  Non-Spherloal  Geometries 

The  result  expressed  in  (3«9)  was  obtained  by  a  Taylor's  series 
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expansion  which  is  valid  for  sufficiently  weak  plasmas  *  Consequently,  since 
the  harmonic  series  representation  of  the  input  admittance  is  absolutely  conver¬ 
gent  for  any  antenna  size,  even  though  it  is  useful  for  computation  only  for 
antennas  which  ara  not  large  relative  to  the  wavelength,  it  can  be  used  to  demon¬ 
strate  the  validity  of  (3, .9)  for  spherical  antennas  of  any  size. 

For  other  antenna  shapes  for  which  the  antenna  contour  is  a  coordinate  sur¬ 
face  in  a  separable  coordinate  system  ( i„e*»  one  in  which  the  wave  equation  is 
separable) ,  a  similar  type  of  representation  of  the  Input  admittance  in  a  harmonic 
series  can  be  carried  out*  It  seems  rather  certain  that  in  such  cases  a  result 
similar  to  (3*9)  would  be  obtained  by  the  same  kind  of  procedure  as  that  followed 
above o 

Re-entry  vehicles  generally  have  rotational  symmetry,  but  seldom  is  the  shape 
''separable't-,  One  is  led  to  conjecture  whether  a  result  similar  to  (3°9)  would 
apply  in  such  cases*  This  would  seem  to  be  worth  investigating  further  by  means 
of  a  generalized  analysis  We  may  express  the  feeling  that  a  proportionality  to 
V  for  small  \>  still  will  be  obtained* 

The  importance  of  a  result  of  the  type  (3*9)  for  upper  atmosphere  probing 
is  worth  pointing  out  in  some  detail.  One  of  the  objects  of  analyzing  the  effects 
of  a  plasma  sheath  on  antenna  impedance  is  to  be  able  to  use  in-flight  impedance 
measurements  to  deduce  the  properties  of  the  plasma*  This  has  led  to  attempts  to 
simulate  the  plasma  so  that  pre-flight  calibrations  can  be  made  of  the  dependence 
of  the  plasma  on  antenna  impedance*  Since  a  plasma  reduces  the  dielectric  constant 
(or  real  part  of  the  refractive  index)  below  unity,  simulation  of  this  character¬ 
istic  is  not  a  simple  matter* 

For  weak  plasmas,  on  the  other  hand,  (3,9)  shows  that  simulation  in  the  strict 
sense  is  not  necessary*  The  validity  of  (3*9),  in  fact,  is  not  limited  in  any  way 
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to  plasma  sheaths,  for  which  if  has  a  negative  real  part,  but  is  valid  for  any 
sheath  for  which  if,  be  it  positive  or  negative  real,  or  complex,  is  small. 

Hence  one  may  use  a  positive  real  if  to  obtain  a  calibration  measurement.  Now  it 
is  a  simple  matter  to  cover  an  antenna  with  a  sheath  of,  say,  polyfoam,  and  thus 
achieve  a  small  known  value  of  y.  By  measuring  the  resulting  change  AYt  a  deter¬ 
mination  of  K,  or  calibration,  of  the  antenna  system  is  achieved. 

3«2  Interference  Phenomenon 

In  Sec.  2.2.3  it  was  pointed  out  that  the  curves  of  AY  vs.  x*  In  Tig. 
2-4  were  suggestive  of  an  interference  phenomenon  between  waves  reflected  from 
the  inner  and  outer  boundaries  of  the  plasma  sheath.  A  qualitative  analysis  of 
this  phenomenon  will  be  developed  here. 

The  geometry  is  illustrated  in  Pig.  3-1,  where  the  region  between  a  and  b  is 
occupied  by  the  Insulating  layer  of  propagation  constant  k,,  and  the  region  be- 


Fig.  3-1  Geometry  for  Analysis  of  Interference  Phenomenon 

tween  b  and  c  is  the  plasma  of  propagation  constant  kg,  while  to  the  right  of  c 
is  free  space,  represented  by  propagation  constant  Arf. 

Assuming  z-  0,  we  have  from  (2.2) 

kg*ngk%  -  (|-x*)%  »-x  CkM/%»  (3.10) 
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Consequently  the  reflection  coefficient  at  the  boundary  e  of  waves  impinging 


from  the  left  is 


where 


ks~k«  „  1 *■  _ &JLz@ 

£  «  Zf.-nKU*~sf\ 


(3.11) 

(3.12) 


At  x2  =  2 „  $  =  45°,  so  that  2@  =  90°c  Consequently  the  reflected  wave  then  is 
shifted  90°  in  phase  relative  to  'the  incident  wave. 

For  Xs >1,  ks  is  pure  imaginary «  as  shown  by  (3,10)  o  Consequently  the  in¬ 
cident  wave  in  advancing  'to  the  right  from  3g  as  well  as  the  reflected  wave  in 
retreating  to  the  left  from  s,  undergoes  attenuation  without  phase  shift.  The 
attenuation  in  a  one-way  passage  between  b  and  e  is 

At  f  =  14  me;,  c~b  ~  0,25  m}  and  x3  =  2  this  amounts  only  to  0.64  dbo  Consequently 
the  insulating  layer  effectively  is  terminated  by  a  pure  reactance  for  X*  ~  2, 
and  this  affects  the  input  admittance  markedly,, 

It  is  quite  evident  from  (3.11)  and  (3.12)  that  the  critical  value  of  x% 
for  a  planar  configuration,  is  2,  since  the  value  of  X3  controls  the  phase  of 
the  reflected  wave.  For  a  non- aero  value  of  z,  the  phase  of  the  reflection  coef¬ 
ficient  is  affected  (as  can  be  seen  from  (3=2)),  and  also  there  is  an  additional 
phase  shift  and  an  additional  attenuation  in  the  wave  travel  between  4>  and  e„ 

For  a  spherical  geometry,  there  is  a  small  departure  from  the  relations  dis¬ 
cussed  above.  The  curves  of  Fig3r,  2-4  and  2-5  clearly  show  the  development  of 
the  phenomenon  as  a  function  of  xa  and  3. 


3.3  Large  Spheres 

The  harmonic  series  representation  of  the  input  admittance  as  expressed 
in  (2.1)  has  been  shown  in  Sec.  2.2,6  to  be  Impractical  for  computational  purposes 

•p 


for  spheres  whose  circumference  is  a  large  number  of  wavelengths.  In  this  Sec¬ 
tion  a  procedure  will  be  outlined  by  means  of  which  it  appears  possible  to  de¬ 
rive  an  alternative  expression  for  the  input  admittance  which  is  useful  for  large 
spheres.  The  procedure  involves  the  following  steps t 

(a)  Transformation  of  (2.1)  into  an  integral  by  means  of  the 
Watson  transformation; 

(b)  splitting  up  this  integral  into  two  components,  one  of  whioh 
corresponds  to  the  geometric  optics  field,  the  other  to  the  diffraction  field; 

(c)  evaluation  of  these  integrals  by  appropriate  methods. 

The  essence  of  the  Watson  transformation  consists  in  recognizing  that  (2.1) 
is  the  sum  of  the  residues  of  a  complex  integral  taken  over  a  suitable  path  en¬ 
closing  the  poles  of  the  integrand  along  the  positive  real  axis.  This  Integral 
is  easily  seen  to  be 


Y- 


f _ «*** _ 

ttolvatl  6„ij,  I  (y*-j4)c osWjt 


(3.13) 


where 


»  n  + 

Z?<*, o)/zf'(*,.), 

z£l'M/z*Hk,a.h 

The  contour  C  is  shown  in  Fig.  3-2. 

The  quantities  -j.*  and  involve  spherical  functions  of  order  n,  or  the 
ordinary  cylinder  functions  of  order  n-*-&=  V. 

In  order  to  put  (3.13)  in  a  form  more  suitable  for  evaluation,  we  split  the 
integral  into  two  parts  by  writing 
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*  l+xtonv'Tr. 


cosyk 

Than  (3.13)  splits  Into  two  parts,  Yt  and  Yz  » 

y,-x -S-VW, 

in  which 

=  <tn>7i  < 

i»«  J  Pj?£(coi0)*in*9d0, 
**-*/* 

V?  mPJT^LcPsB)  d&. 


The  integral  in  Yj  has  no  poles  within  contour  C .  Consequently,  asymptotic 
values  may  be  inserted  for  the  cylinder  functions  in  and  .  Because  of 
the  branch  point  atU=V  (u  denoting  the  argument)  in  the  cylinder  functions  en¬ 
tering  ^ and  /%*,  the  appropriate  asymptotic  forms  for  these  will  have  to  be 
selected . 

In  developing  the  asymptotic  values,  the  usual  asymptotic  expressions  for 
the  associated  Legendre  functions  generally  are  used.  In  [4]  we  have  shown 
that  these  functions  may  be  replaced  by  a  Bessel  function  series  which  is  rapidly 


(3.14) 

(3.15) 
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convergent  for  large  spheres .  This  representation  appears  to  offer  greater 
accuracy  than  the  asymptotic  expansions  of  the  Legendre  functions  for  spheres 
of  intermediate  size. 

In  the  most  general  case,  Yt  may  be  expected  to  be  decomposable  into  several 
components,  corresponding  to  "lens"  and  "rainbow"  terms  due  to  corresponding 
action  of  the  plasma  sheath.  In  almost  all  practical  cases  arising  in  the  re¬ 
entry  problem,  however,  these  terms  would  not  be  expected  to  be  of  any  great  im¬ 
portance  because  of  attenuation  in  the  sheath. 

The  term  Yz  can  be  evaluated  by  residues,  by  deforming  the  contour  C  into 

mi' 

loops  around  the  poles  of  the  integrand,  which  are  the  zeros  of  Z^z^  LA, a)  and 
Ml 

Zy.£  ik,a)  .  These  residues  correspond  to  the  field  which  has  diffracted  around 
the  sphere  one  or  more  times.  Consequently,  it  is  to  be  expected  that  Yt  will  be 
small  compared  to  Y, . 

We  have  given  only  an  outline  of  what  appears  to  be  a  workable  procedure  for 
the  large  sphere  case,  since  the  results  of  the  computation  program  became  avail¬ 
able  only  toward  the  end  of  the  contract  period.  The  procedure  outlined  follows 
in  a  general  way  the  procedure  [5]  which  has  been  found  to  yield  satisfactory  re¬ 
sults  for  the  computation  of  fields  in  the  optical  and  shadow  regions. 


4.  NON-UNIFORM  PLASMAS 


4-1  Introduction 

In  all  the  preceding  work,  the  plasma  has  been  idealised  by  assuming 
both  its  refractive  index  and  thickness  to  be  constant „  Both  of  these  assumptions 
are  far  from  the  true  situation  in  re-entry  plasmas «  Consequently  it  is  desirable 
to  extend  the  treatment  to  a  more  realistic  approximation .  In  this  connection, 
knowledge  of  the  plasma  shape  and  properties  probably  is  still  only  in  a  crude 
state o  Nevertheless,  a  more  realistic  approximation  than  that  of  a  homogeneous 
plasma  should  not  be  difficult  to  obtain. 

In  this  Section,  the  previous  treatment  will  be  refined  to  take  into  account 
an  inhomogeneous  plasma,  but  with  the  restriction  that  the  stratification  be 
spherically  symmetrical,  so  that  the  thickness  of  the  plasma  is  constant.  In  this 
case  it  turns  out  that  the  previous  results  are  modified  by  the  replacement  of  the 
radial  functions  (spherical  Bessel  functions)  by  functions  appropriate  to  the  vari-* 
ation  of  refractive  index  with  the  radial  distance.  Beyond  this,  an  important 
difference  enters,  in  that  the  radial  functions  for  the  electric  modes  are  not 
the  same  as  the  radial  functions  for  the  magnetic  modes,  since  they  satisfy  dif¬ 
ferent  differential  equations. 

A  rather  general  treatment  of  the  problem  of  propagation  in  a  non-homogen eous 
spherically-stratified  medium  was  given  by  B.  Friedman  [6].  Friedman's  analysis 
was  aimed  at  propagation  through  the  earth's  atmosphere.  Consequently,  he  was 
able  to  neglect  certain  terms  arising  from  the  gradient  of  refractive  index  which 
are  small  in  that  problem.  In  the  re-entry  plasma  problem,  however,  these  terms 
no  longer  are  negligible.  In  addition,  Friedman  treated  the  case  of  point  electric 
or  magnetic  dipole  sources  at  a  finite  distance  above  the  surface.  Here  we  shall 
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be  interested  in  the  distributed  surface  source  representing  an  excited  slot. 

In  view  of  these  differences,  it  seems  desirable  to  formulate  the  analysis 
rather  fully,  in  order  to  bring  out  the  differences  between  the  homogeneous  and 
the  inhomogeneous  plasma.  This  will  be  carried  to  the  point  where  the  results  par¬ 
allel  those  for  the  homogeneous  plasma,  so  that  similar  procedures  then  can  be  fol¬ 
lowed. 


4.2  Formulation 

The  underlying  geometry  will  be  the  same  as  in  the  homogeneous  easel 
a  slotted  perfectly-conducting  sphere  is  energized  by  a  voltage  source  at  Its 
center.  The  sphere  is  covered  with  a  uniform  insulating  layer  of  thickness  b-a. 
and  propagation  constant  k, .  Overlying  this  is  a  spherically  stratified  plasma  of 
varying  refractive  index  Viz  =  -nx(R) ,  of  thickness  C-b°  The  plasma  terminates  in 
free  space,  either  abruptly  if  c  is  finite,  or  gradually  if 

As  in  the  homogeneous  case,  the  fields  may  be  separated  into  electric  and 
magnetic  modes.  The  electric  modes  may  be  derived  from  an  electric  Hertz  vector 
potential  T£e  in  virtue  of  the  Maxwell  equation 

div  |g  =  0 . 

This  Hertz  vector,  however,  does  not  satisfy  the  wave  equation.  But  by  writing 

U*-kzSP*  (4.1) 

then  it  may  be  shown  that  Ps  satisfies  the  modified  wave  equation 

*  O,  (4.2) 

where 

(4.3) 

The  quantity  Tig,  may  be  considered  as  an  effective  refractive  index  for  the  electric 
modes. 
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T 

The  electric  and  magnetic  field  etrex-qths  of  the  electric  modes  are  derivable 
-  from  the  equations 

j  E'B-fecuHcurljr*  =Zm  +  <irad(-fejg-),  (4-4) 

tt C m curl  U*  =  'ifc  curl -  (4»5) 

|  it  should  be  noted  that  in  these  equations  kt  appears  rather  than  ^ . 

The  magnetic  modes  are  due  to  source  distributions  which  are  characterized  by 
1  div  %  =  0.  Circular  current  distributions  have  this  property.  The  fields  of  mag¬ 

netic  modes  are  derivable  from  a  radial  Hertz  vector  potential  jTTH,  which,  like 
JTa>  does  not  satisfy  the  wave  equation.  But  by  writing 

JT*  -  (4-6) 

it  turns  out  that  PH  satisfies  the  ordinary  wave  equation 
!  V2PH  +  tfP**  O.  (4°7) 

|  The  electric  and  magnetic  field  strengths  of  the  magnetic  inodes  are  derivable 

from  the  equations 

/?H  =  ^  curl  JTh  *  curl  {klgP*))  (4°8) 

j  H*  a  [Jt > 4radM?kF)\ 3  “  st?  (4o9) 

^  The  advantage  of  writing  the  Hertz  vectors  7fe,  JJ*  in  terms  of  the  scalar 

|  functions  P*  ,  P*  ,  respectively,  is  that  the  latter  quantities  are  solutions  of 

scalar  wave  equations  whose  separation  properties  are  well  known.  Thus,  (4.2) 

!  may  be  separated  in  spherical  coordinates  by  writing 


1 

PPm  =  T(&)Ue(P)  V(<p)  , 

(4»10a) 

1 

where 

1 

! 

~~  +mxV  **  O, 
d^r 

(4.10b) 

1 

i  M*!nB  ie )  +  (s*  “  JlS?)7 *  °> 

(4.10c) 

34 


7j?r  +  (£  -  fi)  t/e  »  o .  (4.ioa> 

The  boundary  conditions  demand  the  continuity  of  the  #  -  and  ^-components  of 
E  and  H  at  the  boundaries  of  the  plasma.  These  lead  to  the  requirement  that 

.  ..z  I  JiM?) 


and 


kl  JR 


be  continuous  across  the  boundaries » 

Similarly,  for  the  magnetic  modes  one  obtains  separable  solutions  of  (4°7)  by 


writing 

KP*  =  T(6)U*(/?)V(f),  (4.11a) 

where  7  and  V'  are  solutions  of  (4.10c),  and  (4.10b),  respectively,  as  before,  but 
U*(k)  is  a  solution  of 


+  ^  U*  *  O.  (4.11b) 

Thus  it  is  clear  that  the  magnetic  radial  functions  d/w  satisfy  a  somewhat  dif¬ 
ferent  differential  equation  than  the  electric  radial  functions  t/*  The  differ¬ 
ence  resides  in  the  appearance  of  the  modified  radial  propagation  constant  k2 
in  (4„10d),  while,  in  (4.11b)  kx  appears  unmodified. 

For  the  magnetic  modes,  the  boundary  conditions  require  the  continuity  of 

U*  and 

aero3s  the  boundaries  of  the  plasma. 

If  one  chooses  the  separation  constant  5^  to  be 

s*  --  n(«+ 1) 

then  a  finite  solution  of  (4.10c)  is 

r*  P„m(cos&). 

Consequently,  a  harmonic  series  representation  for  the  fields,  and  ultimately  of 
the  input  admittance  Y ,  is  obtained  similar  to  that  in  (2.1) .  The  essential  dif- 
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ference  lies  in  the  radial  .functions  for  medium  (2) .  Instead  of  the  spherical 
Hankel  functions  which  appear  (2.1) ,  the  corresponding  solutions  of  (4.10d)  appear. 

For  the  magnetic  modes,  the  same  type  of  harmonic  series  representation  is 
permissible,  the  only  difference  being  that  the  radial  functions  are  the  appropriate 
solutions  of  (4. lib). 

Thus  the  fields  are  derivable  from  the  radial  Hertz  potentials 

7 T*  =  J  £o  at  At  UfWZrCcosS)  cosni,  (4.12a) 

7TH  -  % kl  U«(R)  pnm (cosff)  smmcf.  (4.12b) 

The  field  components  are  obtained  from  (4-4)  and  (4-5)  for  the  electric  modes  and 
from  (4«8)  and  (4.9)  for  the  magnetic  modes.  The  tangential  components  are 


at 


..(*g-^TI  p" 

r*  Hi 


4*  Q*  U * 
*  unm  un 


cosm<f% 


£<f=~fczz{at  7*r'^  sm  wf, 


(4.13) 


unK-%&-atu»'p/)smmri 


mBT 


4 TP  zz{azmk*u?p/-  aLu*'$P$co,m<p. 

In  the  above,  (k^Un)'  denotes  —{ktUn)sn&  '  denotes  * 

Equations  (4.13)  are  valid  for  the  inhomogeneous  (but  spherically-stratified) 
plasma,  medium  (g) ,  as  well  as  the  homogeneous  insulating  layer,  medium  (l) ,  and 
the  outer  free  space  region  (5) «  Only  the  radial  functions  Un  »  U*  will  differ 
for  each  region.  The  forms  of  the  functions  in  each  region  are  determined  by 
the  differential  equations  which  they  satisfy  and  the  boundary  conditions. 

In  order  to  determine  the  fields  by  means  of  the  boundary  conditions,  we 
assume  the  following  forms  for  Li*  and  U*  ; 


4r  »  /£(*,/?)  -Rnjn(k,R\ 


U*1  =  ti\k,R)  - 


(4.14a) 
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4.tZ=  T*- (vf' -  pnirf'l 


TnH(vr™~e"vrf'), 


(4.14b) 


Uf**'cfh?(k*X\  U*^z*ti*>(ksR).  (4.14c) 

iff*  and  iff'  are,  respectively,  outgoing  and  incoming  wave .  solutions  of 
(4-olOd)  for  medium  (§),  while  wf'  and  tirnl"  are  the  corresponding  solutions  of 
(4«Ub) .  and  are  reflection  coefficients  at  the  boundaries  of  the 

sheath,  and  7j;e>H  and  are  the  corresponding  transmission  coefficients. 


4.3  Determination  of  Amplitude  Coefficients 

E  u 

The  amplitude  coefficients  GL„m  and  d„m  are  determined  from  the  pre¬ 
scribed  boundary  values  of  £g  and  £4  on  the  conducting  sphere  R-  0.1 
Ef(a,0,<p)  =  O, 

Egio.,9,^  *  J-  Z^\Spcospf  [/f(0~£ +£)  ~H[e~ ?-£)] 

where  J/Od  is  the  Heaviside  step  function.  From  the  reciprocity  theorem  and  the 
orthogonality  properties  of  the  trigonometric  and  associated  Legendre  functions 
[7],  these  coefficients  are  found  to  be 


(4.15) 


&nm  ~  ~  2n<n+b"  JWuO  IS cosrrup  d$d<p 

"  -10 

=  g"t0)Ja 


rHn+i)  kiUf'b)  ’ 
~  irdWlf  VFCfi 


1  mp,r 


-iro 


&,f)  c o$mfd#d<p 


—  Vw  m  PfX0)X m 

~  n(rw-t)  U^Ujl)  ~  ' 


where 


f-i 


py(ccss)  ,  „ 


(4.16a) 


(4.16b) 


and 
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-JL 

"  2  S  J  Jgm 


Ceos  6) 


?JT(.0) 


de 


and 


r  I,  w  ®  O, 
l£,  Tfi  +  O* 

Except  for  the  functions  6{,e(a)  and  £/nH£a),  (4.16a)  and  (4.16b)  are  the  same 
as  in  the  case  of  the  homogeneous  plasma. 

By  introducing  the  values  of  (4.15a)  and  (4.15b)  into  (4.13),  expressions 
for  the  field  components  in  terns  of  the  exciting  voltage  applied  to  the  slot  are 
obtained.  These  then  can  be  used  to  solve  for  the  fields  and  the  radiation  and 
Impedance  properties  of  the  antenna. 


4„4  Determination  of  Reflection  and  Transmission  Coefficients 

Application  of  the  boundary  conditions,  i.e.  the  continuity  of 

kUE,  UH  ,  and  U*' 

k* 

at  the  boundaries  of  the  sheath,  leads  to  the  following  values  for  the  reflection 
coefficients  j 


B= 


Pi 

*/  0,| 

Ci  Pt 


|E,  6 

C,  P, 


where 


A  _  kfht  h' 

■'  p-v.«j2 

ht  / 1 


» ,  k?h,  hf  I 

*  fcwa. 

n  -  *■'  J 

* s  iww4 

^  •|»SfU‘%  tW-d 


(4.17a) 

(4.17b) 


12  W.n  J»l 

In  terms  of  these  reflection  coefficients,  the  transmission  coefficients  are 
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**  v»*  ’  ' 

Similarly,  for  the  magnetic  modes 
h  _ 

e--v^ 

(Ass 

X>  **  K‘l  Paj 

£*  F*{ " 

j  Cg.  $a| 

-j-H._  ktUu—Xn*j.) 
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where 


hi  A?  .  j 

b£* 

jlv?%  f,Ac%\ 
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£^4 
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fA('%  £*<% | 

1 

1 

it  J/  .  | 

lKMh  lvU*l%\ 

(4,18a) 

(4ol8b) 

(4 -19a) 

(4»19b) 

v4»20a) 
(4° 20b) 


4«5  Input  Admittance 

The  input  admittance  may  be  formulated  in  the  same  way  as  in 


the  case  of  a  homogeneous  plasma.  The  admittance  is  obtained  from  the 

ten 

^  * 


"f~  Jvlfijji  ff £&(&,&,  <f>)  H9ia.,8,<p)  az  &in&  dd  dq. 


-no 


The  result  is 


y*zz-±  r^f-ai  f-e.ifr.wjJ*.  .»  _i fflffimafl 

■  -  il»lj nlml)  "  nfn+O  J’ 


[2]  for 
relation 

(4*21) 

(4«21a) 


where 


®  •-  faijihsL  ~  jo ) 
U-&J*  * 
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(4»22a) 


(4,22b) 


ht-RX+ 


(4»2l)  is  entirely  analogous  to  (2,1),  the  only  difference  being  in  the  re¬ 
placement  of  the  radial  functions  <<$•*,  ^  (defined  by 


(4o23a) 

(4,23b) 

as  in  Sec,  3-3)  by  the  functions  and  ,  respectively. 


4.6  Radiation  Field 

The  radiation  field  is  obtained  by  replacing  the  radial  functions, 
and  by  their  asymptotic  values „  Since  these  radial  functions  are  the  same 

as  in  the  homogeneous  case,  except  for  the  values  of  the  transmission  coefficients, 
the  form  of  the  radiation  pattern  is  essentially  the  same  in  both  caseB,  Some 
important  differences  can  arise,  however,  from  the  variation  of  the  transmission 
coefficients  with  direction,,  In  same  oases,  cut-off  can  set  in  at  certain  angles 
due  to  substantially  total  internal  reflection  at  the  plasma  sheath ,  TMb  phenom¬ 
enon  can  be  affected  appreciably  by  the  nature  of  the  refractive  index  gradient » 
This  phenomenon  also  can  take  place  to  a  different  degree  for  the  electric  and 
magnetic  modes,  due  to  the  fact  that  the  corresponding  radial  functions  for  these 
modes  satisfy  different  differential  equations ,  The  difference  lies  in  the  appear¬ 
ance  of  the  effective  refractive  index  in  the  radial  equation  for  the  electric 
inodes,  while  the  equation  for  the  magnetic  modes  contains  the  refractive  index 
This  difference  will  be  discussed  in  greater  detail  below. 


4.7  Effective  Refractive  Index  for  Electric  and  Magnetic  Modes 

To  illustrate  the  significance  of  the  difference  between  rf2  and  7?2(, 
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consider  the  case  of  a  plasma  which  has  a  refractive  index  which  is  zero  at  its 
inner  radius  R  =  b ,  which  grows  to  a  maximum  at  a  radius  Rm  and  then  decays  (say 
exponentially)  to  zero.  This  can  be  represented,  say,  by 

n*  =  *-***'*. 

The  maximum  value  of  TXz  occurs  where 


Rw-h  =  \/S>. 

Hence 


or,  on  putting 


nl  =A  U?-b)  ~b\ 

R-b  =■  /*, 


=/  re 


(4.24) 


This  variation  is  illustrated  in  Fig.  4-1  by  the  solid  curve. 
To  find  viz*  we  have 

*rMrK\ 

JZ*  r*44"  2F*l 

Hence 


„  d£(  er^\  _  i  rn  ,  y*  g, 

z  dOe*  A*7/*-  /  ^L(r  /^)  +^tj. 


Since  the  bracketed  quantity  is  always  positive,  71%  is  always  less  than  nf .  This 
is  shown  by  the  dotted  curve  in  Fig,  4-1,  Whereas  71%  is  everywhere  finite,  77% 
is  -»  at  r  =  0,  Thus  =  X/rz*  in  (4.10d)  has  a  singularity  at  /*=  0,  whereas 
kx  -  7i%/£  in  (4.11b)  is  regular  at  r  =  0. 

As  in  the  above  example,  the  magnitude  of  the  difference  between  Wf  and  77% 
depends  on  the  "curvature"  of  the  reciprocal  refractive  index  measured  in  terms 
of  the  free-space  wavelength.  Thus  at  sufficiently  high  frequencies  (small  k3) 
this  difference  becomes  small  except  very  near  r  =  0  (the  inner  boundary  of  the 
plasma),  where  the  model  given  by  (4.24)  has  discontinuous  derivatives. 
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-1.0 


Of  course8  the  above  example  is  an  extreme  one  which  is  not  likely  to  be 
encountered  in  practice..  A  more  plausible  assumption  would  be  a  case  where  ti2~ 1 
was  zero  at  R  =  b  and  A?  =  00 o  Nevertheless,  in  any  model  in  which  ti*  becomes  zero 


at  any  point,  nz  becomes  infinite  there o  This  is  evident  from 

1 


This  is  certainly  infinite  unless  s  Z 

As  a  result  of  the  difference  in  the  differential  equations  satisfied  by  Ps 
and  P*l  it  la  to  be  expected  that  there  will  be  a  significant  difference  in  the 

t  IIM- 

behavior  of  the  electric  and  magnetic  modes  which  compose  the  total  field  when  tt* 
differs  greatly  from  Tig,, 


4°8  Discussion 

In  the  foregoing  portions  of  this  section,  the  general  relations  for  the 
fields  and  input  impedance  have  been  developed  for  an  arbitrary  spherically-strati¬ 
fied  plasma  surrounding  the  spherical  slot  antenna,.  It  has  been  shown  that  a  har¬ 
monic  series  representation  of  the  fields  similar  to  that  obtained  in  the  case  of 
a  homogeneous  plasma  can  be  developed  in  this  case,,  The  difference  between  these 
developments  for  the  inhomogeneous  cases  lies  in  the  form  of  the  radial  functions o 
This  results  in  corresponding  differences  in  the  reflection  and  transmission  coef- 
ficientso  Consequently  such  characteristics  as  the  input  admittance,  radiation 
efficiency,  optimization  procedures,  and  radiation  pattern  (especially  at  large 
angles  to  the  direction  of  principal  radiation)  may  be  expected  to  be  affected  by 
the  inhomogeneity  of  the  plasma,.  The  extent  of  these  effects  will  depend  on  the 
nature  and  magnitude  of  the  ^homogeneity „  It  was  hoped  that  data  on  plasma  in¬ 
homogeneity  would  be  available  for  study  of  these  effects  in  connection  with  the 
present  work,  but,  unfortunately,  this  information  was  not  available  to  us« 
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Some  general  remarks  can  be  made  regarding  the  effects  of  an  inhomogeneous 
plasma.  In  Sec.  3.1  it  was  shown  that  a  weak  homogeneous  plasma  produced  a 
change  in  input  admittance  that  is  directly  proportional  to  V  =  fag-  |).  The 
proof  was  based  on  a  Taylor  series  expansion  of  the  radial  functions  which  in¬ 
volve  the  plasma  characteristics.  In  the  case  of  a  weak  Inhomogeneous  plasma,  a 
similar  result  will  be  obtained  in  those  cases  where  the  corresponding  radial 
functions  (both  for  the  electric  and  magnetic  modes)  can  be  expanded  in  a  rapidly 
converging  Taylor's  series.  Thus  one  can  conclude  that  the  result 

AY  =  JLK^>  +•  <Xv*> 

obtained  in  Sec.  3.1  also  holds  for  sufficiently  weak  inhomogeneous  plasmas.  The 
constant  K,  of  course,  will  depend  on  the  nature  of  the  inhamogeneity,  in  general. 

In  [2]  it  was  shown  that  the  external  efficiency  of  the  antenna  in  the  presence 
of  a  homogeneous  plasma  is  the  product  of  three  factors  t  the  ratio  of  power  expended 
in  the  radiating  modes  to  that  delivered  to  all  the  modes,  the  power  transmission 
coefficient  for  the  radiating  modes,  and  the  transmission  factor  at  the  outer 
surface  of  the  sheath  (representing  reflection  loss  at  the  outer  boundary).  In 
the  case  of  an  inhomogeneous  plasma,  the  first  of  these  factors  should  remain  the 
same,  while  the  other  two  could  be  appreciably  different.  An  appreciable  difference 
would  be  expected  in  cases  where  the  overall  depth  of  the  plasma  was  much  greater 
in  the  inhomogeneous  case.  For  a  given  integrated  conductivity  (i.e.  , 

t  -  thickness),  the  attenuation  would  be  about  the  same  in  both  cases,  but  the 
transmission  factor  in  the  inhomogeneous  case  would  be  increased  appreciably  from 
that  in  the  homogeneous  case  (if  the  latter  were  small)  if  the  plasma  extended  out 
to  about  a  body  radius. 
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5o  VOLTAGE  DISTRIBUTION  ALONG  THE  SLOT 


5ol  Introduction 

In  all  of  the  analysis  of  the  center-fed  slot  antenna,  it  has  teen 
assumed  that  the  voltage  distribution  along  the  slot  is  known „  For  the  case  of 
a  slot  which  is  sufficiently  short  compared  with  the  wavelength,  the  voltage  dis¬ 
tribution  cannot  depart  appreciably  from  a  triangular  shape,  either  in  the  pres¬ 
ence  or  absence  of  the  plasma,.  In  the  numerical  computation  program  discussed 
in  Seco  3,  a  triangular  distribution  was  assumedo 

The  assumption  of  a  triangular  distribution  should  be  a  good  one  for  suf¬ 
ficiently  low  frequencies o  At  higher  frequencies,  however,  it  is  to  be  expected 
that  the  voltage  distribution  will  depart  considerably  from  a  triangular  shape  „ 
Furthermore,  the  distribution  may  change  when  the  plasma  sheath  forms „  Conse¬ 
quently,  a  method  for  determining  the  voltage  distribution  or  for  calculating  the 
input  admittance  without  requiring  the  assumption  of  a  known  voltage  distribution 
would  be  desirable,,  In  this  Section,  the  efforts  toward  this  end  wi.ll  be  described,, 
These  efforts  were  not  entirely  successful,  but  a  procedure  which  appears  to  hold 
promise  is  outlined* 


5  o  2  Formulation  of  the  Problem 

The  exterior  admittance  of  a  center  fed  spherical  slot  antenna  was  de¬ 
veloped  in  [2],  with  the  result 

v‘  -Z  L  ‘If  R''m‘r 

z  r'M 


L*i%”(o)i, hT 


Interchanging  the  order  of  the  summations,  this  may  be  written  in  the  form 


(5.1) 
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y#  _  y  Vm  *y« 

r  "mTiWi  *« 


(5.2) 


where 


^  (5.3) 

In  [3],  a  procedtire  was  given  for  calculating  the  quantitea  Z^Kkoi/Z^fki**-) 
and  Zt!i‘\k,Ci)/z^i‘Lk,a.)  for  values  of  n  up  to  about  30  „  Consequently  we  may  as¬ 
sume  that  approximate  values  of  are  available)  at  least  for  the  smaller  val¬ 
ues  of  n  o 

From  the  expression  for  the  azimuthal  component  of  the  magnetic  field,  it 
is  readily  seen  that  this  may  be  written  in  the  form 

//<f>  —  J^0  Hepm  costruf, 

where  H^m  is  independent  of  <p  <,  Then  it  follows  without  difficulty  from  (4.21) 
that 


Y*- 


_  I 


LV(0)3 


it  if 

1  ff  £#  (a,  0,  if)  (a,  9t  <p)  a1-  s  in  $  J&  d<p 


O 


jivco]7  J (fl/ ^ (a> 6'  f) a* 5//7^  cosnupctedf.  (5.4) 

Introducing  the  value  of  *  (&, Qf)  from  (4*15)  ,  it  follows  on  performing  the 

<p  -integration  that 


ym  =  Hfm  *•»*  M- 

?-* 


(5.5) 


In  the  foregoing,  V*  is  the  trt coefficient  in  the  Fourier  expansion  of  the 
slot  voltage 


V(<p)  a  Z  Vm  co»mq>, 

m^O  ' 


(5.6) 


Up  to  now  only  the  exterior  admittance  of  the  slot  has  been  considered. 
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Presumably  if  a  sufficiently  simple  configuration  is  assumed  in  the  interior  of 
the  sphere,  it  would  be  easy  to  obtain  the  interior  field  and  interior  admittance 
of  the  slot  in  much  the  same  manner  as  was  used  for  the  exterior .  In  this  case, 
the  total  admittance  of  the  slot  would  be 


where 


£♦* 


y~  =  =  jzfc *  £  f  (<m  -  h;„)  e  de. 


(5.7) 


(5.8) 


j 

l-i 


Now  the  mih  Fourier  coefficient  of  the  surface  current  density  at  the 

boundary*  between  the  exterior  field  and  interior  field  is  given  by  the  discon¬ 
tinuity  in  the  mode  of  the  magnetic  field  strength 

Km  (0)  ~  ~  Hflns 

Substituting  the  above  into  (5=8) 

M 

Ypr, c  ^  K„  (0)  sin  0  do , 


(5.9) 


or 


where 


15 


y' 

Tt»  “  tmV„  A*» 


V* 


(5.10) 


(5.11) 


Since  the  slot  is  as&l«Jd  to  be  narrow  (s«a) ,  the  above  integral  is  ap¬ 
proximately  equal  to  the  171th  Fourier  coefficient.  Km  (0)1  ,  of  the  current 

density  at  the  equator. 


*  The  boundary  1b  considered  to  extend  entirely  around  the  equator  of  the 
sphere  even  though  the  slot  extends  only  part  of  the  way  around  it.  The  positive 
direction  of  surface  current  flow  (and  also  of  slot  voltage)  is  taken  to  be  the 
minus  9  direction. 
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(5olO)  relates  the  voltage  across  the  slot 


V(f)  ~  2VW 

to  the  surface  current  density  at  the  equator  of  the  sphere 

K(<f)  -  2  Km  costrxpy  (5.12) 

through  the  known  quantities  Ym<,  Consequently  (5.10),  (5.12),  and  (5.6)  should 
be  sufficient  to  determine  the  input  admittance  completely  without  making  any 
assumptions  regarding  the  voltage  distribution,  since  it  should  be  possible  to 
determine  the  latter  through  these  equations. 


5»3  Derivation  of  an  Integral  Equation  for  the  Current  Density 

An  integral  equation  for  K(q>)  will  be  obtained  in  this  section  through 
use  of  (5»10),  (5.12)  and  (5.6).  If  the  mode  radiation  impedance,  Zm*  is 
defined  to  be  the  reciprocal  of  the  mode  radiation  admittance,  (5.10)  can  be 
written  as 

(5.13) 

Since  Km  is  the  tnth  coefficient  in  the  Fourier  expansion  (5.12)  of  KCcf)  ,  one 
has 

ir 

Km  =  ~  jK(<p ')  cos 'd<p '  (5.14) 

Substituting  (5.14)  into  (5.13)  and  multiplying  both  sides  by  cos  /rtf,  there 
follows 


fT 

Vm  cosmcf  =  Zm  a.  caswfj K(ntf')  cos  nup'dtp'.  (5.15) 


By  summing  (5.15)  from  /n=  0  to  °°,  and  interchanging  the  order  of  summation 

and  integration,  (5«15)  becomes 

ir 


w  =  Atz  m  cos  m<p  cos K(q>') a. (5.16) 
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This  equation  holds  for  all  values  of  <jp  in  the  range  ~rl  f  £,  if.  The  in¬ 
tegral  on  the  right  hand  side  of  the  equation  may  be  split  into  two  parts,  one 
representing  integration  from  -  <p,  to  +  <p  along  the  slot,  while  the  rest  repre¬ 
sents  integration  along  the  conducting  surface  of  the  sphere.  Over  the  face  of 
the  slot  the  current  density  must  be  zero  everywhere  except  at  the  center  where 
it  is  fed.  Here  we  set 


K(f')  =  Ioi  SC<p'\ 

where  S(q>')  is  the  Dirac  delta  function  and  I0  is  the  impressed  current  .  Then 

4>. 

r .  *0  go 

J  [I^Zm  c osnup  cos  adcp’  -  &Z  Zm  ccsmtp 

(5.16)  then  can  be  written  as 


VCf)=Iofzm  costrttp  ¥  2/ll>  Ctttruf  ccf/np']K(<p'/'Zcfy'1  (5.17) 
where  use  has  been  made  of  the  fact  that  the  integrand  in  (5.16)  must  be  an  even 
function  of  <p\  from  symmetry. 

If  the  summation  in  (5«17)  is  written  as 


6(<R<p')  5  ZbZm  ccs  mtp  cas  mf‘  (5.18) 

then  (5.17)  becomes 

«r 

V(ep)  =  la  G (f,  O)  ¥  Za  J 6 (#, <f/)  K (cpf)  dtp (5.19) 

When  <p  lies  in  the  range  <f,  <  l<pt  <  m,  Vt<p)  is  zero  since  the  conducting 
strip  in  back  of  the  slot  (see  Fig.  5-1)  is  assumed  to  be  a  perfect  conductor. 
Consequently  (5.19)  becomes 

0  «  GCftO)  ¥  / d<p'  (5.20) 

ft  ^ 

if  <f,  <  Sf  I  <  ir,  while  if  <f>  =  0 

Zfr  =  =  6(0,0)  +  /  S Co,q>')  dtp'.  (5.21) 
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Equations  (5,20)  and  (5,21)  together  with  (5 ,13)  complete  the  mathematical 
formulation  of  the  problem,  (5c 20)  is  a  Fredholm  integral  equation  of  the  first 
kind  with  a  symmetric  kernel 3  GCfprf' )  „  The  unknown  is  the  quantity  ZaKOp )/Io» 
which  is  Za  times  the  current  density  par  unit  impressed  input  current.  If 
(5=20)  can  be  solved,  the  input  impedance,  Zlit  ,  will  be  given  by  (5,21) ,  The 
kernel,  GC<p,<f>')  ,  [as  cast  be  seen  by  substituting  S(q>'~  f  “J  +  ~  '+  <f“)  for 

Kty')  .in  (5,16)]  is  numerically  equal  to  half  of  the  voltage  at  y)  due  to  two  unit 
source  currents,  symmetrically  placed  at  +  <py  and  It  is  necessary  to  use 

symmetrically  placed  sources  in  (5,16)  since  K((f)  was  taken  to  be  an  even  func¬ 
tion  of  <p  when  it  was  expanded  in  a  cosine  aeries.  If  this  assumption  had  not 
been  made  an  integral  equation  would  have  been  obtained  with  a  kernel 
T, Z*,  c&J  m  fp-qv>\  ,  This  kernel  is  a  Green's  function,  and  6 (cp,  cp')  is  the  even 
part  of  this  Green's  function, 

5,4  Solution  of  the  Integral  Equation  for  Some  Special  Cases 

One  standard  method  for  the  numerical  solution  of  integral  equations 
is  to  approximate  the  integral  by  means  of  a  summation,  and  then  to  solve  the 
resulting  set  of  linear  algebraic  equations.  Physically  this  corresponds,  in  the 
case  treated  here,  to  replacing  the  conducting  strip  shown  in  Fig,  5-1  by  a  number 
of  connecting  wires  as  shown  in  Fig,  5-2,  Tile  conducting  surface  has  been  re¬ 
placed  by  ZN  shorting  wires  located  at  <p  q  =  ±  cpxi  ,  0  ,  <p  ="<Pa/  St- ±  K  ° 

The  current  density  in  the  region  f,  £?c  caa  now  be  written  as 

A t  ,, 

Kty  ~  2L  ~£  £ ,  (5.22) 

where  1/  is  the  current  in  the  wire  at  f  =  <pj* 

5,4=1  Formal  Solution 

Substituting  (5-22)  into  (5?20)  and  (5«2l)»  and  using  the  notation 
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FIG.  5-2 


TOP  VIEW 
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(5.23) 


<Pj)  S  G xj 

with  (fa's  0  one  obtains  for  i-  1,  2,  .  .  .  N 


O  =  +  0t  &i-  , 

^'h  =  ^o<>+  0L  <*>y  -^pp- » 


(5.24) 

(5.25) 


or  In  matrix  form 


“<**,  Gw  *  *  •  Go* 

1 

Z2k 

*io  fin  '  '  Sis 

2I/X, 

O 

V 

• 

— 

• 

• 

• 

Gm  *  *  *  sw*_ 

22; A.. 

0 

(5.26) 


Solving  for  Z,;h  by  multiplying  by  the  inverse  of  the  matrix  of  the  equations 
[  \j  ]  there  results 

jGt,  64,  *  *  •«*"! 


_ |>*  <*,-  *  *  VI 

4r" 


<Sn 

•  •  »  •  “  a 

•  *  f  *1  « 

»  •  e### 

G*i4e’**«*« 


(5.27) 


where 


<S*y  =  m^ccsmff.  (5.28) 

Equations  (5.27)  and  (5.28)  provide  a  formal  solution  for  the  input  impedance 
for  the  physical  situation  depicted  in  Fig.  5-2.  This  solution  is  correct,  how¬ 
ever,  only  if  the  matrix  [ <3rxy  ]  in  (5.26)  has  an  inverse,  i.e.  if  the  numerator 
of  (5.27)  does  not  vanish. 

Let  us  consider  the  case  where  the  element  defined  in  (5.28)  is  approxi¬ 
mated  using  Af+  1  modes  j 
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(5.29) 


<3uy  *  z«  °°*  <*>*>»?;  • 

Then  the  matrix  [<?4y  ]  can  be  written  as  a  matrix  product 


Goo  ‘  *  *  64„ 

)  M*4l 

•  •  •  ootMfi 

Z*G  O  •  *•  0 

*i  i 

» 

Q|0  •  *  •  G  ,M 

I  c 0*4, 

•  •  •  emMfi 

ez,d"<o 

*  * 

= 

1  • 

•  »  #  « 

»  •  • 

•  #  • 

•  * 

«  FI  # 

S  •  • 

•  *  e 

'  *  •  cosMfy 

o  o  •  •  • 

(5.30) 


as  can  be  verified  by  taking  the  indicated  product .  Now  the  first  matrix  is  of 
order  /V+  1  by  M+  1,  If  /V  is  greater  than/ty,  then  the  determinant  of  the  product 
vanishes,  [See  W.  L.  Ferrar,  "Algebra"  Oxford  University  Press,  1957,  Theorem  15]. 


5.4.2  Impedance  of  a  Slot  &tgndjng_SggpleteLy  Around  the  Sphere 

One  special  case  to  which  the  formal  solution  (5.27)  may  readily 
be  applied  is  that  of  a  slot  extending  all  of  the  way  around  the  circumference 
of  the  sphere.  The  conducting  strip  in  back  of  the  slot  is  reduced  to  a  narrow 


atrip  located  at  <ft  -  <r r ,  aa  shown  in  Fig.  5-3.  In  this  case  the  strip  evidently 


is  equivalent  to  a  single  wire  at  <f  =  rr  ,  or  to  a  single  pair  of  symmetrically 


placed  wires  at  <f  *  ±  ir.  Taking  ^sf^ciir  (5.27)  becomes 

|Gwo  Izz*  I 

z-  -  *G>,e  6tl‘  r 

"  IG«d  ZZ*(-0am 

,  zSSz~z»t|-<-<>,"4"3 
£  z„ 

ZJZt  +  £z,+Z.Z9+’^^  +  Z,ZM^■ZlZ4.  +  "■  » 

Zo  +  Z)  ■*‘Za+’ 


(5.31) 


5.4=3  Tar^dano  for  Two  Pairs  of  Wires 

Setae  additional  light  is  shed  on  the  general  problem  of  the  slot 
backed  by  a  conducting  strip  if  the  input  impedance  of  a  slot  backed  by  two  pairs 
of  symmetrically  placed  shorting  wires,  one  pair  located  at  fan  and  the  second 
pair  at  f  -  <P,  t  is  compared  with  the  case  already  calculated  above  in  Sec.  5.2, 
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where  the  slot  was  backed  by  a  single  pair  at  <p  a  it  only.  These  two  solutions 
for  the  input  impedance  should  approach  each  other  as  <p,  approaches  it.  That 
this  is  true  only  if  some  resistance  is  Included  in  the  shorting  wires  will  now 
be  shown. 


If  each  shorting  wire  is  replaced  by  a  resistance  R,  then  VCfj)  is  equal 
to  -RIj ,  and  equations  (5.24  and  (5.25)  become 

-*■%'***  & 


s“*  %  S' 


The  solution  is 


6  oc  Got 

(fin+z)  6i* 

•  •  I  *  • 

G»>o*  *  '  * 


Q*hh+j) 


few*!)  6.x 

Gw  (6^f)‘  ‘ 

* 

*  <3tM 

In  the  case  of  shorting  vires  at  f  =  and  s  <yN  •  ir,  this  becomes 


Goo  6*i  G#* 

<*»o  W*t|)  <a,% 


(Gii-*--J)  G*.  I 


Expanding  as  a  polynomial  in  ft/z 
I  Goo  Gm  g. 


ZtA-n  — 


Gio  6h  Gw 
&m  Gw 


*  *  f 


G»»  Am. 


6««  6|*  R .  p*. 


(5.32) 


(5.33) 


(5.34) 


(5.35) 
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As  f  the  first  terms  is  the  numerator  and  denominator  approach  zero,  since 
each  determinant  o on tains  equal  rows*  If  R  is  small,  hut  not  aero,  the  last 
term  may  be  neglected*  R  cancels  out  of  the  ratio  of  the  middle  two  tens  leav¬ 
ing  an  expression  identical  to  the  first  of  equations  (5*31) .  The  limit  of  the 
indeterminate  form  represented  by  the  ratio  of  the  two  leading  terms  of  the  series, 
however,  can  be  shown  not  to  approach  (5 >31) 

This  suggests  that  the  resistance  of  the  conducting  strip  plays  a  vital  part 
in  the  determination  of  the  input  impedance,  and  that  the  approximation  used  in 
obtaining  (5.27),  namely  that  the  resistance  of  the  conducting  strip  is  aero, 
must  be  made  with  care. 

5.5  Attempted  Solution  of  the  Integral  Equation 

One  method  for  obtaining  a  solution  for  the  input  impedance  ifc,  in 
(5.21)  In  the  general  case  (where  the  conducting  surface  is  not  approximated  by 
wires)  is  to  approximate  the  kernel  in  (5. 20)  and  (5.21)  by  &  degenerate  one, 

<3 (<P,  <?')  -  Z  Zm  costnf 

where  the  simulation  now  extends  over  a  finite  number,  M%  of  the  modes  instead  of 
to  infinity.  This  is  the  natural  method  to  use  since  only  a  finite  number  of 
nodes  can  be  calculated  in  any  case.  When  this  is  done  the  integral  equations 
reduce  to  a  set  of  linear  algebraic  equations. 

If  the  order  of  summation  and  integration  is  reversed  in  (5.20)  and  the 
resulting  equation  is  multiplied  by  cos  rup  and  integrated  from  f,  to  -rr»  one 
obtains 

*f-f  M-/ 

(5.36) 

Using  the  same  procedure,  (5.21)  reduces  to 
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(5.37) 


M-t  rt-i 

'  JtZm  +  ZoZm*'*' 

where  the  quantities  xm  and  Tn»  are  defined  by 

_  f  ZaKOfO  j„, 

Xm  *  J  cosnuf  - Y - 

^‘«r 

Tim*  J  cosmf' cosnq>'d<p*. 

Inspection  of  (5.36)  shows  that  a  solution  is  given  by  xw=  -1  so  that 

-  LZm  -  ZZm  »  O.  (5.38) 

This  result  obviously  is  incorrect.  The  reason  why  an  erroneous  result  is 
obtained  by  the  above  approach  will  be  Bade  clear  in  the  next  section  where  the 
problem  will  be  formulated  in  terns  of  infinite  matrices. 


5.6  Matrix  Formulation 

The  problem  may  be  summarized  by  the  five  equations 


vof)  *  Zv„u„($ 

(5.39) 

(5.40) 

V*  *  ( •  +  ionJZm  a-Km 

(5.41) 

V(<?)  *  O ,  <f,<l<pK*- 

(5.42) 

K(f >)■*  T?f(<0, 

(5.43) 

where  the  orthonormal  functions  Unl<f )*  are  used  for  mathematical  con¬ 

venience. 

(5«4l)  can  be  written  in  matrix  font  as 


Vo 

• 

v. 

O  z,  *  .  . 

K, 

• 

*  a 

0  0  z*  .  . 

e 

• 

e 

• 

•  • 

(5.44) 
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or 

V  =  CL-ZV. 

To  express  (5«42)  and  (5.43)  In  matrix  notation,  use  is  made  of  the  idem- 
potent  matrix  A  which  has  the  following  property!  A  is  an  infinite  symmetric 
matrix  which,  when  it  prenrulti  plies  the  column  matrix  of  the  coefficients  of  the 
Fourier  expansion  of  a  given  function,  produces  a  new  set  of  coefficients  which 
correspond  to  a  new  function  identical  with  the  old  in  the  region  M<?|»  tout 
equal  to  zero  if  <f,<  |<pf  <1T*  In  other  words,  multiplication  by  the  matrix 
annihilates  the  portion  of  the  function  which  lies  on  the  back  of  the  slot. 

The  derivation  of  the  elements  of  A  is  as  follows! 


Let 


hW  *  ,  i<ri< 

*  o, 

Let  the  Fourier  expansions  of  -fA(q>)  and  fGp)  be 
fA(q>)  s 

Hep)  *  £ofw  UmOp). 

Then 

♦ft  t^t 

-  / K<p)  Um(<p)  d<p  «  f  Z  K  UnW  Um(<0  dq>, 
-f.  -4i  " 

ft 

/  ur  Up)  Um(<p)  dq> 

ft 

or  in  matrix  notation 


(5.45) 


where 


♦am 


.-L 

-TT 


cos  m<p  cos  r»g>  • 
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Integrating, 


;  I 

zZdf  ~  <jpi/rr, 

-f. 

a  (*'  I  j  v'T  si*nd>i 

**"  J,  ’T7Tcos>'<fd<t  — rr~i 

Anm  *  J  trxp  costupdep  =-^r£ 


s!n(m+n)<fi  ,  s'in(m-y>)<f>r\ 


m+n 


♦  *] 

If  the  function  ^  (q>)  is  defined  by 

V?)  **<<?)> 

s  O, 


then  it  follows  from 


that 


*  f(<p)  -  fA^), 


ni-n 


so  that 


(5.46a) 

(5.46b) 

(5.46c) 

(5.46d) 


(J-A )f  *6f, 

where  I  is  the  identity  matrix  and  3  is  an  idempotent  matrix  that  annihilates 
that  part  of  a  given  function  that  lies  on  the  slot. 

To  write  (5<>43)  in  matrix  form  one  also  needs  the  Fourier  expansion  of 

*<*)  *j|*m6<m6p) 


Then 

it 

=  !  S((?)Um(q>)  ~Um(0) 

-IT 

Thus 


(5.48) 
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In  matrix  notation,  then,  (5.27),  (5.28),  and  (5.29)  become 

M-aZ'K  \ 

BV-O 


(5.49) 


AK=%-$ 

where  a  and  are  scalars,  V  and  K  are  column  matrices  whose  components  are  the 
coefficients  in  the  expansions  (5.39)  and  (5.40),  £  is  a  column  matrix  with 
components  given  in  (5.48) ,  Z'  is  the  infinite  diagonal  matrix  displayed  in 
(5.44),  the  matrix  A  has  components  given  by  (5.46),  and  the  matrix  &  equals  the 
identity  matrix  minus  A . 


Equations  (5.49)  complete  the  matrix  formulation  of  the  problem. 

The  unknown  matrix  V  can  be  eliminated  from  (5°49)  by  substituting  into 
the  second  equation,  giving 


BZ'K^O  \ 


If  in  (5.50)  we  let  K=  Kg+  IoS/a  we  get 

bz'Kb^-^-bz's  " 
AKg  =  o. 


(5.50) 


(5.51) 


The  matrices  in  (5.50)  are  infinite,  so  it  is  not  strictly  correct  to  Bpeak 
of  their  rank  or  of  matrix  inversion.  Roughly  speaking,  however,  it  is  evident 
that  the  matrices  A  8  are  in  a  sense  singular  matrices,  for  in  (5.50)  there 
are  twice  as  many  equations  as  unknowns,  at  least  if  all  of  the  matrices  involved 
are  truncated  after  Af  inodes,  say. 

It  is  now  interesting  to  write  the  integral  equations  (5.20)  and  (5.21)  in 
terms  of  infinite  matrices i  They  become 
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Hero 


(5.52) 

(5.53) 


0=-^ -3Z'S  +  PZ'K* 

Zln  =  *$Z'S  +  Z£iZK, 


Ketti  *  f  Utr,tq')KOv')dQ'  (5*54) 

<tl 

corresponds  to  the  quantities  xm obtained  in  (5.36) ,  where  the  notation  used  is 


slightly  different.  (5.51)  can  be  written  as 

bz'[kb  +  ■#*]»  o. 

This  leads  to  the  solution 

Kr  —  2“ & 


(5.55) 


(5.56) 


which  corresponds  to  the  solution  xM  =  -1  obtained  earlier.  When  the  infinite 
matrix  BZ'  is  truncated  by  using  only  M  inodes,  as  is  done  when  solving  the  in¬ 
tegral  equation  using  a  degenerate  kernel,  the  truncated  matrix  is  non-singular 
and  (5.56)  is  the  only  solution  to  (5.55).  It  happens,  however,  that  the  trun¬ 
cated  matrix  is  very  nearly  singular,  i.e.  its  determinant  approaches  zero  quite 
rapidly  as  the  number  of  modes  used  is  increased.  This  has  been  verified  numer¬ 
ically  using  determinants  of  order  2,  3*  and  4. 

The  nature  of  the  difficulty  encountered  in  the  solution  of  (5.20)  at¬ 
tempted  in  Sec.  5  is  now  clear.  In  approximating  the  infinite  matrix  BZ' by  a 
finite  one  involving  M  modes,  the  singularity  of  the  matrix  has  been  destroyed 
and  all  of  the  non-trivial  solutions  to  the  homogeneous  set  of  equations  (5.55) 
have  been  lost.  There  is  one  more  point  that  should  be  observed!  the  solution 
(5.56)  is  not  consistent  with  (5.54)*  since  S  is  a  matrix  that  corresponds  to 
a  function  that  vanishes  in  the  region  ept<  \<p\  <TT ,  while  (5.54)  requires  that 
the  coefficients  correspond  to  a  function  that  vanishes  in  the  region 

<  <ft.  The  above  discussion  makes  it  clear  that  the  solution  we  seek  is  a 
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nan-trivial  solution  to  the  homogeneous  equations  (5-55)  that  also  satisfies 
requirement  (5*54) •  But  this  solution  is  precisely  the  solution  to  the  pair  of 
equations  (5«5l)» 

The  procedure  that  should  be  used  to  solve  the  problem  is  now  clear.  The 
required  solution  satisfies  (5»5l).  If  (5.51),  which  contains  infinite  matrices, 
is  approximated  using  M  modes,  it  reduces  to  a  set  of  ZM equations  in  M  unknowns. 
It  is  not  expected  that  these  equations  will  be  exactly  consistent.  They  should 
be  approximately  so,  however,  and  a  numerical  solution  should  be  possible.  One 
method  that  suggests  itself  for  solving  (5.51)  both  numerically  and  theoretically 
is  to  use  the  method  of  least  squares.  The  equations  are  first  weighted  by  mul¬ 
tiplying  by  sill  table  diagonal  matrices  Wg  and  WA  • 


‘W.0Zr 

9  • 

'WgBZ'S 

•  **••*■ 

WaA 

Kb 

1 

N 

«  1  t  «  1  ' 

.  O  . 

By  premultiplying  by  the  transpose  of  the  magrix  of  the  equations  (5.57), 
there  results 

[ZW/0Z+  AW^K[Kg  *  Z'BWfBZ'S.  (5.58) 

(5.58)  above  are  the  least  squares  equations  for  (5.51).  The  matrix  of  the 
equations  in  (5.58)  is  easily  shown  to  be  non-singular,  [see  W.  I..  Ferrar, 
loc.  cit.,  Theorem  14],  so  that  (5.58)  may  be  solved  using  Cramer's  rule.  One 
problem  that  must  be  investigated  is  the  proper  selection  of  the  weighting 
matrices  Wa  and  Wg .  That  some  kind  of  weighting  is  necessary  is  obvious  on 
physical  grounds,  since  the  first  equation  in  (5.51)  is  a  relationship  between 
voltages,  while  the  second  has  the  dimensions  of  current.  Also,  the  results  ob¬ 
tained  in  Sec.  5.4.3  suggest  that  it  must  be  necessary  to  replace  (5.50)  by 
something  like 


62 


I 


bz'k^t^bk  ■) 

AK-fs  f  <5-W 

where  ^  is  the  resistivity  of  the  conducting  strip  per  square,  and  is 

the  total  resistance  of  the  strip,  which  oust  remain  small  as  q>t  -pyr« 

5®7  Sunniary 

A  formal  solution  for  the  impedance  of  a  circumferential  slot  in  a 
sphere  which  is  fed  at  one  point  and  short  circuited  at  a  number  of  other  sym¬ 
metrically  located  points  (as  in  Fig,  5-2)  has  been  given  [equation  (5.27)].  A 
simple  formula  (5«3l)  for  one  special  case  (Fig.  5-3)  has  been  obtained.  A  method 
of  attack  is  suggested  for  the  general  case  of  the  impedance  of  a  slot  in  a  sphere 
(Fig.  5-1) o  This  has  not  been  carried  through  as  yet. 
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6c  THE  EFFECTS  OF  PLASMA  RADIATION  ON  RECEIVER  NOISE 

6.1  Introduction 

One  of  the  problems  included  in  the  investigations  conducted  under  the 
subcontract  is  the  analysis  of  the  effects  of  a  re-entry  plasma  upon  receiving 
antennas  aboard  a  hypersonic  vehicle,  with  particular  emphasis  upon  noise  gen¬ 
erating  processes  occurring  in  the  plasma.  This  problem  will  be  dealt  with  in 
this  Section. 

In  the  radiation  of  electromagnetic  waves  from  a  radio  transmitter  within 
a  re-entry  vehicle,  it  has  been  shown  in  previous  reports  that  the  plasma  sheath 
produces  two  principal  effects.  These  are  (1)  an  attenuation  of  the  transmitted 
wave,  and  (2)  a  lowering  of  the  input  impedance.  The  second  effect  in  most  cases 
would  lead  to  an  additional  loss  above  that  in  (l)  due  to  dissipation  in  the  im¬ 
pedance  matching  device  that  would  be  required  between  the  transmitter  and  the 
antenna,  or  to  increased  reflection  losses  if  no  attempt  were  made  to  restore  an 
impedance  match  in  the  presence  of  the  plasma  sheath. 

From  the  previous  analyses,  it  1b  possible  to  calculate  the  field  strength 
produced  at  a  receiving  point  (on  the  ground,  say)  by  a  known  voltage  impressed 
across  the  antenna  terminals  in  the  vehicle.  Since  the  input  impedance  of  the 
antenna  has  been  determined,  the  received  field  strength  also  can  be  expressed 
in  terms  of  the  power  delivered  to  the  transmitting  antenna  terminals.  By  the 
reciprocity  theorem,  therefore,  the  voltage  received  across  the  antenna  terminals 
in  the  vehicle  from  a  ground  transmitter  can  be  determined. 

Although  the  reciprocity  theorem  allows  the  previous  treatments  of  the  trans¬ 
mission  problem  to  be  extended  to  the  reception  problem,  this  applies  only  to 
the  (desired)  signal.  Reciprocity  is  strictly  a  two-terminal  relationship,  how¬ 
ever,  so  that  it  does  not  consider  the  signal  from  any  other  source  which  may  be 
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active  in  the  received  signal  band.  In  particular,  noise  generated  fay  some 
other  source,  including  noise  generated  by  the  plasma,  is  not  taken  into  account 
in  such  a  treatment. 

In  determining  whether  a  radio  circuit  is  capable  of  providing  reliable 
communication  of  intelligence,  a  criterion  is  adopted  for  a  minimum  signal-to- 
noise  ratio.  The  noise,  in  many  cases,  is  that  generated  in  the  receiver  by 
thermal  and  other  fluctuations  in  its  circuits  and  components.  The  magnitude  of 
this  internal  noise  is  dependent  on  the  temperature  of  the  noise-generating  com¬ 
ponents,  which  usually  is  taken  to  be  that  of  the  surroundings.  In  the  case  of 
a  low-noise  receiver,  the  limiting  noise  is  that  received  from  external  sources 
which  radiate  to  the  antenna.  Thus  a  natural  question  arises  when  a  high-temper¬ 
ature  plasma  surrounds  the  antenna  of  the  vehicle!  Is  the  effective  noise  temper¬ 
ature  of  the  receiver  equal  to  that  of  the  plasma? 

In  this  Section  this  question  will  be  examined.  It  will  be  shown,  in  fact, 
that  the  effective  noise  temperature  of  the  receiving  system  depends  not  alone 
on  the  temperature  of  the  plasma,  but  also  on  the  attenuation  through  the  plasma 
at  the  frequency  in  question.  As  a  consequence,  a  lower  optimum  frequency  usually 
will  exist  for  reception  through  such  a  plasma  than  for  transmission. 

6.2  Equivalent  Noise  Temperature 

In  dealing  with  the  noise  contributions  from  various  sources,  it  is 
convenient  to  use  the  concept  of  equivalent  noise  temperature.  It  then  becomes 
a  relatively  simple  matter  to  determine  the  noise  received  by  the  antenna  from 
the  plasma. 

The  situation  is  Identical  to  that  encountered  in  the  field  of  radio  astronomy 
in  the  reception  of  emissions  from  radio  noise  sources.  The  derivation  of  the 
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well-known  relations  will  be  given  here,  both  for  completeness  and  to  bring  out 
the  principles  Involved. 

An  antenna  receives  an  amount  of  power  from  a  radio  aource  which  depends  on 
the  emission  characteristics  of  the  Bource,  the  attenuation  in  the  space  between 
source  and  receiver,  and  the  receiving  characteristics  of  the  antenna.  The  amount 
of  power  received  usually  is  expressed  in  terms  of  an  equivalent  source  tempera¬ 
ture  7*.  This  temperature  is  that  of  a  black  body  which  would  produce  the  same 
received  power  at  the  frequency  of  observation,  From  Planck's  law  of  radia¬ 
tion,  unit  area  of  such  a  black  body  would  radiate  an  amount  of  power  per  unit 

frequency  band  per  unit  solid  angle  equal  to 

.  Zhf  I 

<6.i) 

where  h  is  Planck's  constant,  and  4k.  is  Boltzmann's  constant.  This  radiation  is 
randomly  polarized.  Since  in  most  radio  work  h i/kT«  1,  the  above  equation  re¬ 
duces  to  the  Rayleigh-Jeans  approximation 

E-  (6.2) 

Thus  the  emitted  power  density £"  may  be  expressed  in  terms  of  the  equivalent 
temperature  Tby  means  of  (6.1),  or  the  simplified  fora  (6.2)  if  applicable.  £  is 
also  called  the  brightness  of  the  source . 

The  above  concept  of  an  equivalent  brightness  and  equivalent  temperature 
also  may  be  applied  to  a  coherent  or  modulated  type  of  signal. 

In  observations  of  discrete  sources,  the  source  will  subtend  a  solid  angle 
jQ.  at  the  receiver  which  is  smaller  than  the  solid  angle  of  the  antenna  beam. 
The  received  power  then  is  given  by 

=JtT3-Q-/£U  (6.3) 

where  B  is  the  receiver  bandwidth,  A  =&3/4ris  the  effective  receiving  area  of 
the  antenna,  <5  its  gain,  and  JXA  =  4*/<5.  The  factor  arises  from  the  fact  that 
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the  antenna  responds  only  to  a  fixed  polarization,  whereas  the  emitted  radiation 
is  assumed  to  be  randomly  polarized. 

If  the  signal  received  from  the  desired  souroe  is  to  be  detectable,  it  must 
be  sufficiently  large  to  produce  an  observable  increase  in  the  noise  output  of 
the  receiver  when  the  antenna  is  pointed  at  the  source.  Receiver  internal  noise 
frequently  is  expressed  in  terms  of  an  equivalent  noise  temperature  7#.  The 
noise  output  of  the  receiver  in  the  absence  of  any  incoming  radiation  from  space 
then  Is  the  same  as  if  a  noise  power 


Ph  *  *TnB 


(6.4) 


were  supplied  to  the  input. 

This  noise  power  also  is  frequently  expressed  in  terms  of  the  receiver  "noise 
figure"  F°  This  is  the  ratio  of  the  actual  receiver  output  when  supplied  with 
an  available  power  P0  from  a  resistor  (or  other  noise  source),  at  a  reference 
equivalent  temperature  To,  to  the  output  if  TH  were  zero  (i.e.,  if  the  receiver 
generated  no  internal  noise) .  Since  noise  powers  are  additive, 

p  —  P^'t'  fp  =  7v  t  7p 
Pc  ”  To  ' 


From  this. 


r*=  (f-\)To. 


(6.5) 


The  reference  temperature  T0  usually  is  taken  to  be  290°  K.  F  frequently  is  ex¬ 
pressed  in  decibels. 

It  is  evident  from  the  above  discussion  that  the  addition  of  noiBe  powers 
is  equivalent  to  the  addition  of  the  corresponding  equivalent  noise  temperatures. 

(6.4)  may  be  used  to  define  an  equivalent  temperature  for  the  radiation  re¬ 
ceived  by  the  antenna.  Denoting  this  by  7J,  it  follows  from  (6.3)  that 

Ta  -  T-O/Jla-  (6.6) 

In  the  above  discussion,  it  has  been  assumed  implicitly  that  there  is  free- 
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space  transmission  from  source  to  receiving  antenna.  This  means,  among  other 
things,  that  no  attenuating  regions  between  source  and  receiver  are  supposed  to 
be  present.  The  effect  of  attenuation  will  now  be  considered. 

6.3  Effect  of  Attenuation  on  Receiver  Noise  Temperature 

In  considering  the  effect  of  radiation  from  the  plasma  on  the  noise 
at  the  input  of  the  receiver,  it  is  convenient  to  think  of  the  plaama  as  a  trans¬ 
ducer.  A  transducer  composed  wholly  of  reactance  elements  generates  no  noise, 
since  it  possesses  no  mechanism  for  imparting  random  motions  to  the  electrons 
circulating  in  its  elements.  Noise  generation  in  a  transducer  thus  is  synonymous 
with  dissipation.  A  transducer  having  a  very  large  (-*«>)  loss  would  impart  noise 
of  temperature  T0  to  the  receiver,  where  T0  1b  the  transducer  (ambient)  temperature. 

Consider,  for  example,  a  transmission  line  which  connects  an  antenna  to  a 
receiver.  Let  the  line  loss  factor  be  L  (i.e.,  Pc,*t/Pin  ~ L) ,  and  the  noise  temper¬ 
ature  of  the  receiver  be  In  -  The  transmission  line,  through  its  dissipation,  re¬ 
duces  the  temperature  of  the  signal  passed  to  the  receiver  to  and  also  con¬ 

tributes  noise  of  equivalent  temperature  (l-L)  TL  ,  where  7i  is  the  line  temperature. 
Hence  the  total  noise  input  to  the  receiver  is 

+  0-47L.  <6»7> 

The  signal-to-noise  ratio  thus  is 

(4/4)~4Z/fc+(i-t)ZLJ,  (6.8) 

whereas  in  the  absence  of  line  loss  it  would  be 

(n/p*)i*  Tcl/z, . 

The  reduction  in  signal-to-noise  ratio  thus  is 

l'*  n/r„i<L,  <6-’> 

and  thus  is  greater  than  the  line  loss  factor  itself. 
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Some  numerical  examples  at  this  point  may  be  helpful  in  assessing  the  order 
of  magnitude  of  the  above  effect .  Consider  a  line  having  an  effective  noise 
temperature  of  4000°,  and  a  receiver  with  a  noise  figure  of  6  db  (F  *-4).  Then  from 
(6.5) 

-  3TE  =  870°  K. 

Suppose  the  line  attenuation  at  the  operating  frequency  is  10  db,  i.e.  L  -  0.1. 
Then  from  (6.7)  the  effective  receiver  temperature  is 

Tj  -  870  +  0.9*4000  =  4470°  K. 

Consequently  the  signal-to -noise  ratio  (in  the  absence  of  external  noise)  which 
existed  in  the  absence  of  line  attenuation  is  reduced  by  the  factor  L\  which  from 
(6.9)  is 

L  =  0.1/[1  +  0.9*4000/870]  =  0.0194  — -17.1  db 
On  the  other  hand,  if  the  line  attenuation  is  only  3  db,  then 

Tf!  ~  870  +  0.5*4000  =  2870°  K 

and 

L'  =  0.5/[l  +  0.5*4000/870]  =  0.152  -*--8.2  db. 

In  a  similar  way,  a  plasma  surrounding  a  receiving  antenna  contributes  noise 
to  the  receiver  input.  This  noise  depends  both  on  the  equivalent  temperature  of 
the  plasma  (75 )  and  on  its  attenuation  (/,)  at  the  frequency  to  which  the  antenna 
is  tuned.  Consequently,  in  a  determination  of  an  optimum  frequency  for  reception 
through  a  re-entry  plasma,  the  following  factors  are  involved 1 

(a)  the  signal  transfer  characteristic,  including 

(b)  the  plasma  attenuation; 

(c)  the  equivalent  temperature  of  external  noise  sources; 

(d)  the  plasma  temperature. 
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6.4  Ontlmi  nation  Procedures 


In  optimizing  the  power  radiated  from  a  vehicle  through  a  re-entry 
plasma,  only  factors  (a)  and  (b)  above  are  involved.  For  the  case  of  a  strong 
uniform  plasma,  for  example,  we  have  shown  [2]  that  the  optimum  frequency  is  the 
one  for  which  the  plasma  thickness  is  2J-  skin  depths. 

In  the  receiving  case,  the  transfer  characteristic  for  the  signal  and  for 
external  noise  is  the  same.  Consequently  the  ratio  of  signal  to  external  noise 
is  not  affected  by  the  plasma.  But  in  the  absence  of  external  noise,  the  signal- 
to-ncise  ratio  is  decreased  when  an  attenuating  plasma  is  present,  as  shown  by 
(6.3).  Hence  the  decrease  in  signal -to-noise  ratio  due  to  an  attenuating  plasma 
is  not  as  great  when  external  noise  is  present  than  when  no  external  noise  exists. 

The  level  of  external  noise  is  a  composite  of  contributions  from  terrestrial 
and  cosmic  sources.  Approximate  values  as  a  function  of  frequency  for  various 
locations  and  times  are  available  [8]  for  use  in  planning  or  design  purposes.  In 
view  of  the  great  variability  in  external  noiBe  level,  it  is  not  possible  to  make 
an  exact  specification  of  optimum  frequency  for  all  situations.  However,  a  typi¬ 
cal  procedure  that  can  be  followed  will  be  illustrated  below. 

To  illustrate  the  optimization  procedure,  consider  the  case  of  a  strong  uni¬ 
form  plasma  of  constant  thickness.  For  this  situation,  the  analysis  of  [2]  may 
be  used  to  determine  the  radiated  fields,  etc.  The  various  factors  listed  at  the 
end  of  Sec.  6.3  will  be  considered  in  turn. 

(a)  Signal  transfer  characteristic 

For  the  case  of  a  strong  plasma  sheath,  the  radiation  character¬ 
istics  (including  sheath  attenuation)  are  such  that  the  high-frequency  range  is 
of  interest.  A  slotted  Bphere  antenna  then  has  a  radiation  pattern  equivalent 
to  that  of  a  loop  whose  plane  is  perpendicular  to  the  center  of  the  slot.  The 
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(6.10) 


far  field  then  is  given  by  [2,  p.  22,  eq«  (27a)] 

3ft  c  Vi  e-1**1'  ,  -:k~ic.-b) 

where 

•>(  =  k^Lb-o). 

For  the  usual  case  of  a  slot  length  small  relative  to  the  wavelength,  the  voltage 
distribution  along  the  slot  is  approximately  triangular,  so  that 

The  last  factor  in  (6,10)  represents  the  sheath  attenuation.  Denoting  this 
/  ife 

attenuation  by  L  , 

we  have 

JL 


(6.11) 

(6.12) 


l/r  I  ^  3  e  -g  v(0)  sm&'/'k 

Frl  2- CL  l  r  L  * 

(6.12)  gives  the  free-space  field.  The  actual  field  at  the  receiver  (which  we 

assume  ground-based)  will  be  modified  by  propagation  conditions.  In  order  to  sim¬ 
plify  the  treatment  here,  we  will  assume  free-space  conditions  to  exist  in  the  en¬ 
suing  discuB3ion. 

Let  the  receiving  antenna  have  an  effective  gain  (i.e.,  including  line  losses) 
of  6t,  so  that  its  effective  aperture  is 

A  =  Qrrf+K.  (6.13) 

Furthermore,  assume  the  antenna  is  matched  to  a  50-ohm  line,  so  that  the  load 


m 


Gr 

30 


impedance  is  50  ohms.  Then  the  power  delivered  to  the  load  is 

l5y>|aA  '<* 

^  I207T 

and  consequently  the  load  voltage  is 

^  =  I7^Tv(0>(Sr£)/‘ 

C\5?*cJL 

87 ra  r 

where  the  magnitude  of  K  is  independent  of  frequency.  Hence 
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V,  _  K(6rL) 


tlx 


(6.14) 


m  it^*i  * 

(6«14)  gives  the  voltage  ratio  between  the  slot  and  the  input  to  the  ground- 
based  receiver. 

The  coupling  between  the  antennas  in  transmission  from  the  slot  antenna  to 
the  ground-based  antenna  can  be  represented  as  a  mutual  impedance,  so  that  vre  can 
draw  the  equivalent  circuit  shown  in  Fig.  6-1  (a).  The  mesh  equations  are 

V{0)~  X?>Zn  +  l2'ZM*Xt'>Z*, 


+  irZal  —  V2 


r  <*>  ss _ T< 

Iat  Zzi  +  Zp 


«> 


For  2,,.=  50  ohms, 

Vi.  =  ~  £  WO) 

so  that 

»  ZVzZ»/V(0). 

For  the  reverse  case  of  transmission  from  the  ground-based  antenna  to  a  re¬ 
ceiver  connected  to  the  slot  antenna,  we  can  draw  the  equivalent  circuit  shown 
in  Fig.  6-1  (b) .  In  an  analogous  way  we  obtain 


(6.15) 


3 


Z,» 


Zo 


(6.16) 


Zn Z»  Z»+ZRl‘ 

By  the  reciprocity  theorem,  £,=  J&*  that  from  (6.15)  and  (6.16) 

1/  -  -  r«>2_  *-£-  Zji.  ,_Zg,  _ 

^  ZniZff,  WO)  Zx%  2,,+Zai  • 

Assuming  that  the  receiver  input  is  tuned  optimally  (hy  an  automatic  tuning  device, 
for  example)  so  that 


7* 


the  above  equation  reduces  to 

=  Zn.  Z;! 

v'  ¥10)  Z,»  ZK,  WO)  ioo>e;,|Y(*  Wo)  100#  * 


I 


1 


(6.17) 


where  &  =  <fi>(Y)  is  the  input  conductance  of  the  slot.  In  the  above  we  have  taken 
Z„=  50. 
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(6.18) 


Vr  may  be  expressed  in  terns  of  the  transmitted  power  Pj,  since 

Pr  *  VT*/Z2X. 

Since  the  signal  power  a.%  the  receiver  input  is  we  have 

Ps  ®  (  4j^r|l+AlA)  ^  *  ^L>  (6.19) 

where  ff  may  be  called  the  signal  transfer  characteristic.  is  equivalent  to  a 

signal  temperature  Ts\ 

rs'=  Ps'/AP. 


(6.20) 

(6.21) 


In  the  absence  of  the  sheath,  the  signal  power  would  be 
r  m  PrGr  3*l*‘**8'^  StGt 
™  4-&1  8icrx  *  44F7 A/  *  i 

where  / 

3 

A' ,  8ir  \  r  J  * 

The  frequency-dependent  terms  in  ,5  of  (6,19)  are  6^^  and  j(.  For  a  ground- 
based  antenna  of  constant  effective  aperture,  6,*,/^,  while  J0rmay  be  found 

from  calculated  values  of  Y,  Consequently  we  may  calculate  the  signal  transfer 
characteristic  &  and  plot  a  curve  of  it  as  a  function  of  frequency. 

(b)  Plasma  attenuation 

The  plasma  attenuation  is  given  hy 

4-  |e~*,‘e'w|. 

But  for  a  highly  conducting  plasma 

4  = 

where  <tyis  the  radian  plasma  frequency,  the  collision  frequency,  and  vthe  veloc¬ 
ity  of  light.  Hence 

l  =  e-z&tO-lc-b)  _  e--^(^r)  Cc~^  (6,22) 

/.  thus  decreases  exponentially  as  the  square-root  of  frequency.  Frosts  known  (or 
estimated)  values  of  ajy  and  >/,  L  may  be  calculated  and  plotted  against  frequency. 


(c)  Equivalent  temperature  of  external  noise 

External  noise  is  a  combination  of  terrestrial  (natural  +  man-made) 
and  cosmic  noise.  The  noise  level  is  a  complicated  function  of  location,  time 
(of  day  and  season),  as  well  as  of  frequency.  For  planning  purposes,  curves  pub¬ 
lished  by  the  CCIR  [8]  may  be  used.  These  curves  actually  give  the  equivalent 
noise  temperature,  FA ,  expressed  in  db  relative  to  Tc~  288°  as  a  reference.  Hence 
a  curve  of  external  noise  tanqaerature  vs.  frequency  may  be  selected  for  the  loca¬ 
tion  and  time  of  interest.  This  noise  temperature,  however,  applies  for  a  re¬ 
ceiving  antenna  which  is  a  short  (relative  to  the  wavelength)  grounded  vertical 
rod.  Since  such  an  antenna  has  a  gain  of  3/2  relative  to  an  isotropic  antenna, 
the  CCIR  noise  temperatures  should  be  multiplied  by  2/3  (i.e.,  decreased  by 
roughly  2  db). 

In  making  use  of  the  CCIR  curves,  it  is  necessary  to  assume  that  the  external 
noise  is  uniformly  distributed  over  all  directions.  This  actually  1b  not  true, 
since,  for  example,  natural  terrestrial  noise  is  due  to  thunderstorms,  which  are 
not  uniformly  distributed;  nor  are  propagation  characteristics,  which  are  super¬ 
imposed  on  the  source  distribution,  independent  of  direction.  In  accordance  with 
the  above  assumption  of  a  uniform  distribution  of  external  noise,  the  solid  angle 
XX  of  the  source  is  4ft>  Hence  from  (6.5)  and  (6.6)  the  external  noise  temperature 
at  the  antenna  is 

Te  =  (/=*-!)  To  -  I  (A  -  0  -288  *  432(^a-|),  (6.23) 

where  the  gain  of  the  antenna  has  been  taken  to  be  3/2. 

Since  external  noise  experiences  the  same  reduction  in  passing  through  the 
plasma  sheath  as  does  the  signal,  the  temperature  of  external  noise  at  the  receiver 
may  be  obtained  by  multiplying  (6.23)  by  the  ratio  Ft/Fs  obtained  from  (6.19)  and 

(6,21) i 
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./  -  K&  L 


Te  ~  p,Te  *  '432^“^  (6*24) 

Values  of  7J/  nay  be  calculated  from  (6,24)  and  plotted  as  a  function  of 
frequency. 


(d)  Equivalent  plasma  temperature 

The  equivalent  plasma  temperature  already  has  been  discussed  in 
Sec,  6»3.  If  we  denote  the  plasma  temperature  by  TP,  then  from  (6,7)  the  total 
noise  temperature  of  the  receiver  (excluding  external  noise)  is 

T"  *  rN  +  U-L)TP.  (6,25) 

Since  L  (and  possibly  T^)  is  a  function  of  frequency,  values  of  T/f  may  be  plotted 
vs,  frequency  from  (6,22)  and  data  on 


From  the  plotted  curves  of  the  various  factors  discussed  in  (a)  -  (d)  above, 
it  is  a  simple  matter  to  combine  these  to  obtain  a  curve  of  the  overall  signal- 
to-noise  ratio  at  the  receiver,,  This  is  given  by 


where  is  given  by  (6,19),  by  (6,24),  and  L  by  (6.22),  From  such  a  curve, 
the  optimum  frequency  can  be  determined. 

In  general,  because  of  the  large  effect  of  plasma  attenuation  both  in  re¬ 
ducing  the  signal  level  and  in  raising  the  internal  noise  level,  it  is  to  be  ex¬ 
pected  that  the  optimum  frequency  for  reception  will  be  lower  than  the  optimum 
frequency  for  transmission.  The  extent  to  which  one  can  gain  in  reception  by 
lowering  the  frequency  is  determined  largely  by  the  external  noise  level.  The 
latter  la  a  strong  function  of  geographical  location  and  temporal  factors,  so 
that  the  optimum  frequency  for  reception  will  vary  accordingly. 
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7.0  CCMCLUSIQNS  AND  HBCQMMEMPATIQNS 


7.1  Conclusions 

7.1.1  Calculation  at  Twr^rt. 

A  study  of  calculations  of  input  admittance  of  a  spherical  slot 
antenna  made  by  MSVD  on  the  IBM  704  revealed  two  significant  phenomena  and  two 

limitations : 

(a)  A  linear  relation  exists  between  the  change  in  input  ad¬ 
mittance  fro*  the  froe-space  value  and  the  ohange  in  refractive  index,  as  ex¬ 
pressed  in  the  relation 

for  snail  (7i4-  I ) . 

(b)  An  apparent  interference  phenomenon  takes  place  for  small 
collision  frequencies  around  a  certain  value  of  plasma  frequency. 

(c)  The  calculations  are  erroneous  for  certain  values  and 
ranges  of  the  parameters. 

(d)  The  series  formulation  used  is  not  practical  for  values 
of  vehicle  circumference  greater  than  about  one  wavelength. 


7.1.2  TiwrfttPtl 

The  following  results  wen  obtained  in  extension  of  previous 

developmental 

(a)  The  relation  for  AY  given  above  in  Sec.  7.1.1  (a)  is 
shown  to  exist  for  any  vehicle  else.  This  suggests  a  useful  technique  for  pre- 
flight  calibration  of  the  antenna  ao  that  in-flight  measurements  may  be  used  to 
determine  plasma  properties. 
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(b)  The  phenomenon  mentioned  in  See.  7.1.1.(b)  is  explained 
as  the  interference  produced  by  reflection*!  free  the  outer  boundary  of  the 
pi— »  shsath,  whereby  the  antenna  is  effectively  terninated  by  an  almost  pure 
reactance. 

(o)  An  analysis  of  the  input  adnittsnoe  which  is  suitable  for 
large  vehicle  sizes  is  outlined. 

7.1.3  fttfCTPgflPtPVf  flMBt 

A  general  formulation  is  given  of  the  problem  of  an  inhomogeneous 
spherical  plasma.  It  is  shown  that  the  previous  forms  of  the  result  are  still 
retained,  the  only  change  being  that  the  radial  functions  require  alteration. 

In  this  connection,  a  significant  new  feature  is  that  the  electric  and  magnetic 
modes  satisfy  different  differential  equations . 

7.1.4  Slot  Voltage  Distribution 

The  determination  of  the  voltage  distribution  along  the  slot  in 
the  general  case  leads  to  an  Integral  equation.  A  practical  solution  to  this 
problem  was  not  found,  but  a  method  of  attack  which  nay  prove  fruitful  is  pre¬ 
sented. 


7.1.5  Plasma  Woise 

In  a  discussion  of  the  effect  of  noise  generated  in  the  plasms 
upon  the  reception  problem,  it  is  shown  that  the  effective  noise  temperature 
of  a  receiving  system  aboard  a  re-entry  vehicle  depends  on  the  attenuation  of 
the  plasma  aa  well  as  on  the  plasma  temperature.  Consequently,  the  optimum 
frequency  for  reception  usually  will  be  significantly  lower  than  for  transmission. 
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7.2  Rg£ggag&aaUgg.e. 


1.  Is  view  of  the  inaccuracies  of  the  double  precision  program  used 
In  the  calculations,  it  is  reccsnended  that  a  triple  precision  progran  be  em¬ 
ployed  in  any  future  calculations. 

2.  The  practical  importance  of  a  calibration  technique  based  on  the 
relation  AV  *  a.K(t\v-\)  makes  it  desirable  to  extend  the  validity  of  this  ex¬ 
pression  to  geometries  other  than  spherical. 

3.  In  Seo.  3.3,  a  procedure  was  outlined  whereby  practically  useful 
formulas  can  be  obtained  for  large  spheres.  The  details  of  this  procedure 
should  be  worked  out  so  that  formulas  suitable  for  numerical  computation  will 
be  available. 

4.  In  Sec.  4,  formulations  were  derived  for  Inhomogeneous  plasmas. 

It  is  recosmended  that  these  formulations  be  applied  to  available  information  on 
plasma  properties,  so  that  the  quantitative  effeota  of  plasma  inhamogenelties  on 
antenna  radiation  properties  can  be  deduced. 

5.  The  effect  of  the  plasma  on  the  voltage  or  current  distribution 
along  an  antenna  is  still  an  unsolved  problem,  la  view  of  its  importance  in 
upper  atmospheric  research,  further  work  on  this  problem  is  justified. 

6.  The  effect  of  noise  generated  by  the  plasma  on  reception  aboard  a 
re-entry  vehicle  has  been  shown  in  Sec.  6  to  lead  to  lower  optimum  working  fre¬ 
quencies  in  reoeption  than  in  transmission.  It  is  recommended  that  experiments 
aimed  at  verifying  this  conclusion  be  considered. 
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